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Extended quantum systems can be theoretically described in terms of the Schwinger-Keldysh functional inte¬ 
gral formalism, whose action conveniently describes both dynamical and static properties. We show here that in 
thermal equilibrium, defined by the validity of fluctuation-dissipation relations, the action of a quantum system 
is invariant under a certain symmetry transformation and thus it is distinguished from generic systems. In turn, 
the fluctuation-dissipation relations can be derived as the Ward-Takahashi identities associated with this sym¬ 
metry. Accordingly, the latter provides an efficient test for the onset of thermodynamic equilibrium and it makes 
checking the validity of fluctuation-dissipation relations unnecessary. In the classical limit, this symmetry re¬ 
duces to the well-known one which characterizes equilibrium in the stochastic dynamics of classical systems 
coupled to thermal baths, described by Langevin equations. 
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I. INTRODUCTION 

In recent years, the question under which conditions and 
how a quantum many-body system thermalizes has received 
ever-growing attention. This interest has been primarily trig¬ 
gered by the increasing ability to prepare and manipulate such 
systems, which might be either isolated [1-3] — as it is typi¬ 
cally the case in experiments with cold atoms [4, 5] — or in 
contact with an environment (open), and therefore subject to 
losses and driving. 

After an abrupt perturbation, isolated systems are generi- 
cally expected to thermalize in the sense that expectation val¬ 
ues of local quantities at long times can be determined on the 
basis of suitable statistical ensembles [1, 3]. However, this 
might not be the case because of the presence of an exten¬ 
sive amount of conserved quantities induced by integrability 
[6-10] or because of a breaking of ergodicity due to the oc¬ 
currence of many-body localization [11-14]. Although it is 
possible to define a variety of effective temperatures based on 
the static [15, 16] and dynamic properties [17, 18] under such 
circumstances, the lack of thermal behavior is witnessed by 
the fact that these temperatures do not necessarily take all the 
same thermodynamic value. 

Examples of open systems include exciton-polaritons in 
semiconductor heterostructures [19, 20], arrays of microcavi¬ 
ties [21, 22], trapped ions [23], as well as optomechanical se¬ 
tups [24]. In general it is unclear, a priori, by which physical 
mechanism an effective temperature is possibly established in 
these systems and, in case, what determines its value. Recent 
work, however, suggests possible mechanisms where an ef¬ 
fective temperature can occur as a consequence of the compe¬ 
tition between driven-dissipative and coherent dynamics [25- 
33]. Irrespective of its cause, effective thermalization often 
affects only the low-energy degrees of freedom [16, 25-37]. 

All these examples show clearly that the presence of effec¬ 
tive thermodynamic equilibrium (which might be established 
only in a subsystem or within a specific range of frequencies) 


in the steady state of a system is often by no means obvi¬ 
ous. Hence, before addressing the question of whether the 
time evolution of a certain system leads to thermalization or 
not, it is imperative to identify criteria which allow a clear- 
cut detection of thermodynamic equilibrium conditions in the 
stationary state. In this regard, it is important to consider not 
only the static properties of the density matrix of the system, 
which describes its stationary state, but also the dynamics of 
fluctuations; being encoded, e.g., in two-time correlation and 
response functions, it might or might not be compatible with 
equilibrium. As a fundamental difference between static and 
dynamic properties, the latter necessarily involves the genera¬ 
tor of time evolution, while the former does not. 

In this work we consider the following operative definition 
of thermal equilibrium: a system is in thermal equilibrium at a 
certain temperature T if expectation values of arbitrary prod¬ 
ucts of operators, evaluated at different times, are connected 
by quantum fluctuation-dissipation relations (FDRs) involv¬ 
ing the temperature T. These FDRs were shown [38-40] to be 
equivalent to a combination of the quantum-mechanical time- 
reversal transformation [41] and the Kubo-Martin-Schwinger 
(KMS) condition [42, 43]. Heuristically, the latter condition 
expresses the fact that the Hamiltonian ruling the time evolu¬ 
tion of a system is the same as that one determining the den¬ 
sity matrix of the canonical ensemble, which characterizes the 
system when it is weakly coupled to a thermal bath. In both 
the generalized FDRs and the KMS condition the temperature 
appears as a parameter. 

From the theoretical point of view, static and dynamical 
properties of statistical systems (both classical and quantum) 
are often conveniently studied in terms of dynamical function¬ 
als, which are used in order to generate expectation values of 
physical observables in the form of functional integrals over 
a suitable set of fields. Then, it is natural to address the issue 
of the possible equilibrium character of the stationary state by 
investigating the properties of the corresponding dynamical 
functional. In the case of classical statistical systems evolving 


2 


under the effect of an external stochastic noise of thermal ori¬ 
gin, this issue has been discussed to a certain level of detail in 
the past [44-49], and it was found that the dynamical func¬ 
tional acquires a specific symmetry in thermodynamic equi¬ 
librium. As in the case of the FDRs and the KMS condition, 
the (inverse) temperature j3 - IjT enters as a parameter in this 
symmetry transformation. Remarkably, classical FDRs can 
be derived as a consequence of this symmetry. For quantum 
systems, instead, we are not aware of any analogous deriva¬ 
tion based on the symmetries of the corresponding dynamical 
functional, which takes the form of a Schwinger-Keldysh ac¬ 
tion (see, e.g.. Refs. 50-57). 

The aim of the present work is to fill in this gap by show¬ 
ing that also the Schwinger-Keldysh dynamical functional of 
a quantum system in thermal equilibrium is characterized by 
a specific symmetry, i.e., it is invariant under a certain trans¬ 
formation Tp. This symmetry may be considered as the gen¬ 
eralization of the classical one mentioned above, to which it 
reduces in a suitable classical limit [58]. In addition, Tp can 
be written as a composition of the quantum-mechanical time 
reversal expressed within the Schwinger-Keldysh formalism 
— reflecting a property of the generator of dynamics — and 
of the transformation which implements the KMS conditions, 
associated with a property of the state in question. The ex¬ 
istence of this symmetry was already noticed in Ref. 58 for 
mesoscopic quantum devices, where it was used to derive fluc¬ 
tuation relations for particle transport across them. However, 
to our knowledge, the connection between this symmetry and 
the presence of equilibrium conditions has not yet been estab¬ 
lished. 

The rest of the presentation is organized as follows: the 
key results of this work are anticipated and summarized in 
Sec. II; in Sec. Ill, we specify the symmetry transformation 
Tp, provide its various representations, and list a number of 
properties which are then detailed in Sec. IV. In particular, we 
discuss the invariance of unitary time evolution in Sec. IV A, 
while in Sec. IV B we consider possible dissipative terms 
which are invariant under Tp. We discuss how the quantum 
symmetry reduces in the limit /i —> 0 to the one known in clas¬ 
sical stochastic systems in Sec. IV C. As we discuss in Sec. V, 
the symmetry can be interpreted as a practical implementation 
of the KMS condition on the Schwinger-Keldysh functional 
integral. Finally, Sec. VI presents applications of the equilib¬ 
rium symmetry: in Sec. VIA we derive the FDR for two-point 
functions while in Sec. VIB we show that the steady states of 
a quantum master equation explicitly violate the symmetry. 
The case of a system driven out of equilibrium by a coupling 
with two baths at different temperature and chemical potential 
is considered in Sec. VIC; Sec. VID briefly touches upon a 
number of other applications of the symmetry. 


II. KEY RESULTS 

a. The invariance under Tp of the Schwinger-Keldysh ac¬ 
tion is a sufficient and necessary condition for a system to be 
in thermal equilibrium. As mentioned in Sec. I, we consider 
a system to be in thermal equilibrium if all the FDRs are satis¬ 


fied with the same temperature T - (5^^ or, equivalently [38- 
40], if the KMS condition (combined with time reversal) is 
satisfied. In Sec. V, we show that these conditions imply the 
thermal symmetry 7]g of the Schwinger-Keldysh action cor¬ 
responding to the stationary state of the system. Conversely, 
the fluctuation-dissipation relations can be derived as conse¬ 
quences of the symmetry, proving their equivalence. 

b. A different perspective: thermal equilibrium as a sym¬ 
metry. A key conceptual step forward we take in this work 
is to provide a compact formulation of thermal equilibrium 
conditions of a quantum system — i.e., the KMS condition 
(or, alternatively, of the equivalent hierarchy of FDRs) — 
in terms of a single symmetry Tf, which can be considered 
as the fundamental property of quantum systems in thermal 
equilibrium. This perspective is especially fruitful within the 
field-theoretical formalism, where various tools have been de¬ 
veloped to work out the consequences of the symmetries of 
the action of a given system. In this context, for example, 
the hierarchy of generalized quantum FDRs can be derived 
straightforwardly as the Ward-Takahashi identities associated 
with the thermal symmetry (see Secs. V and VIA). In ad¬ 
dition, the Schwinger-Keldysh formalism provides a conve¬ 
nient framework to take advantage of very powerful and effi¬ 
cient renormalization-group techniques for studying the possi¬ 
ble emergence of collective behaviors and for monitoring how 
the effective description of a statistical system depends on the 
length and time scale at which it is analyzed. The possible 
scale dependence of the restoration/violation of the equilib¬ 
rium symmetry could shed light on the mechanism underlying 
the thermalization of extended systems. 

As we mentioned above, the idea of viewing thermal equi¬ 
librium as a symmetry is certainly not new. However, while 
previous studies were primarily concerned with classical sta¬ 
tistical physics [44-49], here we generalize this idea to the 
quantum case. 

c. Unification of the quantum and classical description of 
equilibrium systems. As pointed out in Ref. 58, the equilib¬ 
rium symmetry reduces, in the classical limit, to a known sym¬ 
metry which characterizes thermal equilibrium in open clas¬ 
sical systems [44^9]. In Sec. IV C we review the classical 
limit of the Schwinger-Keldysh action for a system coupled 
to a thermal bath [55, 56] and we discuss in detail how the 
classical equilibrium symmetry is recovered. The compari¬ 
son with the classical symmetry highlights some remarkable 
differences with the quantum case: in fact, in classical sys¬ 
tems, thermal equilibrium can be regarded as a consequence 
of detailed balance, which, in turn, is related to the property of 
microreversibility of the underlying microscopic dynamics. In 
fact, the classical equilibrium symmetry is derived by requir¬ 
ing the dynamical functional to satisfy these properties [46- 
48]. For quantum system, instead, an analogous satisfactory 
definition of detailed balance and microreversibility is seem¬ 
ingly still missing, leaving open the important question about 
the very nature of thermal equilibrium of quantum systems. 

d. Efficient check for the presence of thermodynamic equi¬ 
librium conditions. The symmetry is of great practical value, 
as it reduces answering the question about the possible pres¬ 
ence of thermodynamic equilibrium to verifying a symmetry 
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of the Schwinger-Keldysh action instead of having to check 
explicitly the validity of all FDRs. In particular, we show in 
Sec. VIB that the Markovian quantum dynamics described 
by a Lindblad master equation [59, 60] explicitly violates the 
symmetry. This reflects the driven nature of the system: in¬ 
deed, the Lindblad equation may be viewed as resulting from 
the coarse graining of the evolution of an underlying time- 
dependent system-bath Hamiltonian, with a time dependence 
dictated by coherent external driving fields. 

Moreover, in Sec. VIC we consider a bosonic mode cou¬ 
pled to two baths at different temperatures and chemical po¬ 
tentials: in this case, the resulting net fluxes of energy or par¬ 
ticles drive the system out of equilibrium with a consequent 
violation of the symmetry. 

e. A new perspective on the construction of the 
Schwinger-Keldysh action. At the conceptual level, the exis¬ 
tence of the symmetry provides a new perspective on the con¬ 
struction of Schwinger-Keldysh functional integrals. In par¬ 
ticular, as customary in quantum held theories, one may con¬ 
sider the symmetry as the fundamental principle: indeed, it is 
explicitly present for any time-independent (time-translation 
invariant) Hamiltonian which generates the dynamics of a sys¬ 
tem at the microscopic scale. Then, requiring the symmetry 
to hold for the full effective Keldysh action at a different scale 
Axes the admissible dissipative terms so as to satisfy FDRs 
between response and correlation functions of arbitrary or¬ 
der; translating back into the operator language, this provides 
a concrete hint why stationary density matrices of the form 
p ~ e^^^ are favored over arbitrary functions p{H) for the 
description of static correlation functions. 


III. SYMMETRY TRANSFORMATION 

As we anticipated above, a convenient framework for the 
theoretical description of the time evolution of interacting 
quantum many-body systems is provided by the Schwinger- 
Keldysh functional integral formalism [55, 56]. It offers 
full flexibility in describing both non-equilibrium dynamics 
and equilibrium as well as non-equilibrium stationary states, 
which is out of reach of the finite-temperature Matsubara tech¬ 
nique [61]. In addition, it is amenable to the well-established 
toolbox of quantum held theory. The simplest way to illus¬ 
trate the basic ingredients of the Schwinger-Keldysh formal¬ 
ism is to consider the functional integral representation of the 
so-called Schwinger-Keldysh partition function Z. For a sys¬ 
tem with unitary dynamics generated by the Hamiltonian H 
and initialized in a state described by a density matrix po, 
this function is given by Z = \i{e~’^'poe'^'^. (Note that, as 
it stands, Z = 1; however, it is instructive to focus on its struc¬ 
ture independently of its actual value.) In this expression, time 
evolution can be interpreted as occurring along a closed path: 
starting in the state described by po, the exponential e^'^‘ to 
the left of Po corresponds to a “forward” evolution up to the 
time f, while the exponential e‘^‘ to its right corresponds to an 
evolution going “backward” in time. The trace tr(- ■ ■) con¬ 
nects, at time f, the forward with the backward branch of the 
time path and therefore it produces a closed-time-path inte¬ 


gral. Along each of these two branches, the temporal evo¬ 
lution can be represented in a standard way as a functional 
integral of an exponential weight e‘^ over suitably introduced 
(generally complex) integration variables, i.e., fields, ifr+{t,x) 
and i/^_(f, x) on the forward and backward branches, respec¬ 
tively. These fields are associated with the two sets of coher¬ 
ent states introduced as resolutions of the identity in-between 
two consecutive infinitesimal time evolutions in the Trotter 
decomposition of the unitary temporal evolution along the two 
branches [55, 56]. The resulting Schwinger-Keldysh action S 
is a functional of if±(t, x) and it is generally obtained as a tem¬ 
poral integral along the close path in time of a Lagrangian 
density. (Explicit forms of S will be discussed further below, 
but they are not relevant for the present discussion.) By in¬ 
troducing different (time-dependent) sources for the fields 

on the two branches, the partition function Z[7+, 7_] is no 
longer identically equal to 1 and its functional derivatives can 
be used in order to generate various time-dependent correla¬ 
tion functions (see, e.g.. Refs. 55-57). 

As we show further below in Sec. V, a system is in ther¬ 
modynamic equilibrium at a temperature T - 1//3, if the 
corresponding Schwinger-Keldysh action is invariant under a 
certain transformation 7^ which acts on the fields along the 
closed time path. In order to specify the form of 7^, we fo¬ 
cus on the dynamics of a single complex bosonic held, which 
is sufficiently simple but general enough to illustrate conve¬ 
niently all the basic ideas. In this case, the transformation Tp 
turns out to be composed of a complex conjugation [62] of the 
field components if/a- with cr = +, an inversion of the sign of 
the time variable, and a translation of the time variable into 
the complex plane by an amount icry6/2, i.e., 

TfHl/a-{t,x) = if^li-t + /cr/?/2,x), 

Tpil/l.(t,x) = il/c-i-t + io-pi2,x). 

For convenience and future reference we provide an alter¬ 
native compact representation of the action of Tjj both in 
the time and real space domain (f, x) as well as in the 
frequency-momentum domain (at, q). The convention for the 
Fourier transforms of the fields, conveniently collected into 
two spinors 'I'o.(f, x) = {\fa-{t,x),\f/*^{t,x)Y, is the following: 

In this relation, d is the spatial dimensionality of the system, 
and the field spinors in the frequency-momentum domain are 
deflned as 'Td-(w, q) = q), iA^(-<u, -q))^. Accordingly, 

we can write the symmetry transformation Tjj in the form 

7>'T^(f,x) = + iaPM, 

7]3'I'^(cu,q) = q), 

(3) 

where we introduced the Pauli matrix cr^ — (i o)- The trans¬ 
formation in real time requires evaluating the fields for com¬ 
plex values of the time argument, which in principle is not de¬ 
flned; however, the complementary representation in Fourier 
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space indicates how this can be done in practice: in frequency backward branches: 

space, the shift of time by an imaginary part icrfill amounts ^ ^ 

to a multiplication by a prefactor = — (i/r+ + i/^_), ^ - (/r_). (4) 

V2 V2 


As usual within the Schwinger-Keldysh formalism, it is 
convenient to introduce what are known as classical and quan¬ 
tum fields. These are defined as the symmetric and antisym¬ 
metric combinations, respectively, of fields on the forward and 


Combining these fields into spinors q) = 

((^v(w,q),0*(-w,-q))^ — where the index v - c,q dis¬ 
tinguishes classical and quantum fields — the transformation 
Tp takes the following form, which we report here for future 
reference: 


T;i<i>c(a),q) = cr^(cosh(J3cL>/2)(^c(-ca,q) - sinh08w/2)cl>^(-w,q)) , 
q) = o-^ (- smh(J3oj/2)^ci-ca, q) + cosh06w/2)O^(-w, q)) . 


We anticipate and summarize here a number of properties of 
the equilibrium transformation 7^, which are going to be dis¬ 
cussed in detail in Secs. IV and V: 

1. The transformation is linear, discrete and involutive, 
i.e., 7^ = 1. The last property follows straightfor¬ 
wardly from Eqs. (1) or (3). Concerning linearity, note 
in particular that the complex conjugation in Eq. (1) 
affects only the field variables, i.e., 7jg/l(/'o-(f, x) = 

-H icrpi2,x) for T e C (see Sec. V B). 

2. 7^ can be written as a composition Tj} - T o 'Kp of a 
time-reversal transformation T and an additional trans¬ 
formation 'Kp, which we will identify in Sec. V C as 
the implementation of the KMS condition within the 
Schwinger-Keldysh functional integral formalism. 

3. Tp is not uniquely defined, due to a certain freedom in 
implementing the time-reversal transformation within 
the Schwinger-Keldysh functional integral formalism, 
as discussed in Sec. V B. However, without loss of gen¬ 
erality, we stick to the definition provided by Eq. (1) 
and we comment on the alternative forms in Sec. V B. 

4. The transformation Tp leaves the functional measure in¬ 
variant, i.e., the absolute value of the Jacobian determi¬ 
nant associated with Tp is equal to one, as discussed in 
Sec. V D and shown in App. D. 

5. The various forms of the transformation Tp presented 
above apply to the case of a system of bosons with van¬ 
ishing chemical potential /u. In the presence of fi ^ 0, 
Eq. (1) becomes 

Tp,^ipAt, x) = H- io-f3l2, x), 

+ /cr/J/2,x), 

with a consequent modification of Eq. (3), which can be 
easily worked out. After a transformation to the basis of 
classical and quantum fields according to Eq. (4), this 
modification amounts to shifting the frequency a> in the 
arguments of the hyperbolic functions in Eq. (5), i.e., to 

W —> W - jU. 


6. In taking the Eourier transforms in Eqs. (3) and (5) one 
implicitly assumes that the initial state of the system 
was prepared at time t - -oa, while its evolution ex¬ 
tends to t - oa. In the following we will work under 
this assumption, commenting briefly on the role of an 
initial condition imposed at a finite time in Sec. IV C. 


IV. INVARIANCE OF THE SCHWINGER-KELDYSH 
ACTION 

As we demonstrate further below in Sec. V, a system is in 
thermodynamic equilibrium if its Schwinger-Keldysh action 
S is invariant under the transformation Tp, i.e., 

5m=5[7j3'P], (7) 

where, for convenience of notation, 'T = (^+, i/^*, 
collects all the fields introduced in the previous section into a 
single vector. The tilde in S indicates that all the parameters 
in S which are related to external fields have to be replaced 
by their corresponding time-reversed values (e.g., the signs of 
magnetic fields have to be inverted), while in the absence of 
these fields the tilde may be dropped. 

According to the construction of the Schwinger-Keldysh 
functional integral outlined at the beginning of the previous 
section, the action corresponding to the unitary dynamics of 
a closed system is completely determined by its Hamilto¬ 
nian H. The initial state po of the dynamics enters the func¬ 
tional integral as a boundary condition: if the system was 
prepared in the state po at the time t - 0, the matrix el¬ 
ement (t/^+,olpo|iA-,o), where |iA±,o) are coherent states, deter¬ 
mines the (complex) weight of held configurations at the ini¬ 
tial time with i^±{Q,x) = iA±,o(x). In Sec. IV A, we demon¬ 
strate the invariance of the Schwinger-Keldysh action asso¬ 
ciated with a time-independent Hamiltonian dynamics under 
the transformation Tp. In particular, this invariance holds 
for for any value of f. Interestingly enough, the Schwinger- 
Keldysh action associated with a Hamiltonian of a simple non¬ 
interacting system — which can be diagonalized in terms of 
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single-particle states — turns out to be invariant under an en¬ 
hanced version of this transformation, involving possibly dif¬ 
ferent values of p for each of the single-particle states (see 
Sec. IV A 3). A constraint on the value of p, however, comes 
from the inclusion of the boundary condition for the func¬ 
tional integral which specifies the initial state po- Here we 
are interested in the stationary state of the system, which is 
generically reached a long time after its preparation in the 
state po- Hence, we assume that this was done in the distant 
past, i.e., at f = -oo, and that the evolution of the system ex¬ 
tends to f = -Hoo (cf point 6 in Sec. III). In the construction 
of the Schwinger-Keldysh functional integral for a system in 
thermodynamic equilibrium [55, 56], a convenient alternative 
approach for specifying the appropriate boundary conditions 
corresponding to the initial equilibrium state po of the sys¬ 
tem, consists in adding infinitesimal dissipative contributions 
to the action. Usually [55, 56], the form of these contribu¬ 
tions is determined by the requirement that the Green’s func¬ 
tions of the system are thermal with a specific temperature 
T - lip, i.e., that they obey a fluctuation-dissipation relation; 
once these terms are included, any reference to po may be 
omitted. We demonstrate in Sec. IV B 1, that these dissipative 
contributions are invariant under 7p with exactly the same p. 
Hence, the thermal symmetry provides a different perspective 
on the construction of the Schwinger-Keldysh functional in¬ 
tegral for a system in thermal equilibrium; while the unitary 
contributions are fixed by the Hamiltonian of the system, the 
requirement of invariance under the symmetry transformation 
7p can be taken as the fundamental principle for specifying 
the structure of the dissipative terms which can occur in the 
action if the system is in thermodynamic equilibrium at tem¬ 
perature T - lip. We emphasize that only the simultaneous 
presence in the Schwinger-Keldysh action of both the Hamil¬ 
tonian and the dissipative contributions yields a well-defined 
functional integral; the dissipative terms in the microscopic 
action are taken to be infinitesimally small as for an isolated 
system, where they merely act as a regularization which ren¬ 
ders the functional integral finite and ensures that the bare re¬ 
sponse and correlation functions satisfy a FDR; on the other 
hand, if the isolated system is composed of a small subsystem 
of interest and a remainder which can be considered as a bath, 
then finite dissipative contributions emerge in the Schwinger- 
Keldysh action of the subsystem after the bath has been inte¬ 
grated out. This scenario is considered in Secs. IV C and VIB. 
Moreover, the system can act as its own bath; in fact, one 
expects the effective action for the low-frequency and long- 
wavelength dynamics of the system to contain dissipative con¬ 
tributions which are due to the coupling to high-frequency 
fluctuations. In Sec. IV B, we explicitly construct dissipative 
terms which comply with the thermal symmetry 7p. In partic¬ 
ular, we And that the noise components associated with these 
dissipative terms must necessarily have the form of the equi¬ 
librium Bose-Einstein distribution function, as appropriate for 
the bosonic fields which we are presently focussing on. 


A. Invariance of Hamiltonian dynamics 

The Schwinger-Keldysh action S associated with the dy¬ 
namics generated by a time-independent Hamiltonian H can 
be written as the sum of a “dynamical” and a “Hamiltonian” 
part, 5dyn and respectively, 

S - Sdyn + S-Ji, (8) 

5dyn = ^ r Kitr,5,'I'+-'I'licr,5,(9) 

^ Jf,X 

= -( 10 ) 

where we used the shorthand f = dt, f = f d'^x, while 
cr^ — (0 ) is the Pauli matrix. This structure of the 

Schwinger-Keldysh action results from the construction of the 
functional integral outlined at the beginning of Sec. III. In 
particular, the Hamiltonians 'N± are matrix elements of the 
Hamiltonian operator N in the basis of coherent states 
i.e., 'Tfo- = (^/ro-|Ff|(Ao-)/(iAcrl<A<T), where the amplitudes of 
the coherent states are just the integration variables in the 
functional integral [55, 56]. Henceforth we focus on the case 
of a bosonic many-body system with contact interaction, i.e., 
with Hamiltonians in Eq. (10) given by 

Here m is the mass of bosons, t the chemical potential, and 
A parametrizes the strength of the s-wave two-body interac¬ 
tion. We consider this case because it is sufficiently general 
for the purpose of illustrating all basic concepts associated 
with the thermal symmetry and, in addition, in the classical 
limit it allows a direct comparison with classical stochastic 
models [63, 64], where <pc - (iA+ + <A-) / V2 plays the role of 
a bosonic order parameter held. This point is elaborated in 
Sec. IV C. 

Below we show that the invariance of the Schwinger- 
Keldysh action S under Tp is intimately related to the struc¬ 
ture of the action, i.e., to the fact that it can be written as the 
sum of two terms containing, separately, only fields on the 
forward and backward branches. 


1. Dynamical term 

To begin with, we show that the dynamical contribution 
S dyn to the Schwinger-Keldysh action S given in Eq. (9), is 
invariant under T^g, i.e., that SdyrXTpdp - 5dyn[‘l’]. To this 
end, it is convenient to express the original fields (iA±,i/'![] 
in the so-called Keldysh basis, which is formed by the 
classical and quantum components {(pc,q,4>*q] introduced in 
Eq. (4). Eor the sake of brevity, we ari'ange these fields into 

the vector ■ Rewriting 5dyn in these 

terms and in frequency-momentum space, we obtain (J^ ^ s 
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f dwd‘'q/( 2 ;r)‘'+') 

5dyn[7>d>] = I w [cosh^(y8w/2)<l)J(w, q)crj<Dc(w, q) 

- sinh^(/?w/ 2 )<l)^(w, q)cr 2 <l)^(w, q) + sinh(; 6 w/ 2 ) cosh(ySw/ 2 ) 
X (Oj(w, q)fr,d)J (w, q) - OjCod, q)%(co, q))]. (12) 

The combination Ov(w, q)cr 20 y(w, q) = 0*(w, q)0v(w, q) - 
(f>vi-co, -q)(f>l(-ai, -q) with v = c, g' is an odd function of 
(w, q), whereas a)smh(J3(jj/2)cosh(J3a>/2) is even, and there¬ 
fore the integral over the product of these terms vanishes. 
Then, with some simple algebraic manipulation, the first two 
terms in Eq. (12) are recognised to be nothing but 5dyn[®], 
from which the invariance of S dyn follows straightforwardly. 
Note that this property holds independently of the value of the 
parameter p in the transformation Tp. 


2. Hamiltonian contribution 

We consider now the transformation of the Hamiltonian 
contribution S-ft in Eq. (10) under T)?. Eirst, we argue that 
the strictly local terms (i.e., those which do not involve spatial 
derivatives) in the Hamiltonian (11) are invariant under 7 ) 3 ; 
then, we extend the argument to the case of quasilocal terms 
such as the kinetic energy contribution oc |Vi/r±p or non-local 
interactions. Consider a contribution to 5./r of the form 

W]= r (v+(f,x)-v.(t,x)), (13) 

Jt,\ 

where Vo-(f, x) = (iA^(f, x)il/a-(t, x))^ is a generic local contribu¬ 
tion to the Hamiltonian TYo- and N is an integer. In particular, 
for At = 1 we obtain the term proportional to the chemical po¬ 
tential in Eq. (11), while for N - 2, W'['I'] is just the contact 
interaction. Since Vo-(f, x) is real, under the transformation Tjg 
[see Eq. (6)] only its time argument is shifted according to 
TpVa-it, x) = Vo-(-f -H icrpi2, x) and, taking the Eourier trans¬ 
form with respect to time of this relation, one eventually finds 

7>v<y(w,x) = e-"'^"/V(-w,x). (14) 

Accordingly, the vertex (13) is invariant under Tp\ in fact, be¬ 
ing local in time, its diagrammatic representation — where 
the fields x) and i/'[j-(f, x) are represented by ingoing and 
outgoing lines, respectively — satisfies frequency conserva¬ 
tion for in- and outgoing lines, as can be seen by taking the 
Eourier transform of each of the fields in Vo-(f, x) individually. 
In particular, the frequency w in Eq. (14) corresponds to the 
difference between the sums of the in- and outgoing frequen¬ 
cies and only the w = 0 component contributes to Eq. (13). 
(As stated above, we assume that the time integrals in Eqs. (9), 
(10), and therefore (14) extend over all possible real values, 
i.e., we focus on the stationary state of the dynamics.) This 
component, however, is invariant under Tp as follows directly 
from Eq. (14), and hence the same is true for the vertex, for 
which ^[7>T'] 


Clearly, the invariance of the vertex and of the dynamical 
term in Eq. (9) relies on the fact that vertices, which are lo¬ 
cal in time, obey frequency conservation. (Note that, as in 
Sec. IV A 1, this invariance holds independently of the value 
of the parameter p in Tp.) Accordingly, one concludes that 
any contribution to the Hamiltonian, which is local in time 
and does not explicitly depend on time, is invariant. In par¬ 
ticular, the proof of invariance presented here for the ver¬ 
tex in Eq. (13) can be straightforwardly extended to expres¬ 
sions containing spatial derivatives such as the kinetic energy 
oc in Eq. (11) and even to interactions which are not 

local in space, as long as they are local in time, as anticipated 
above. Note, however, that these considerations do not rule 
out the possible emergence upon renormalization or coarse- 
graining of terms which are non-local in time, as long as they 
are invariant under Tp. This case is discussed further below in 
Sec. IV B. 


3. Enhanced symmetry for non-interacting systems 

The equilibrium transformation Tp presented in Sec. Ill in¬ 
volves a single parameter p. While this form is appropriate 
for the Gibbs ensemble describing the thermal equilibrium 
state of the interacting many-body system with the Hamil¬ 
tonian in Eq. (11), an enhanced version of the symmetry is 
realized in non-interacting systems. Since these systems can 
be diagonalized in terms of single-particle states, they are triv¬ 
ially integrable. Statistically, integrable systems are described 
by a generalized Gibbs ensemble [ 6 , 8 , 65-71], constructed 
from the extensive number of conserved quantities (with pos¬ 
sible exceptions, see, e.g.. Refs. 67-70). In the case of non¬ 
interacting systems which we consider here (or, more gener¬ 
ally, for any system that can be mapped to a non-interacting 
one), these integrals of motion are just the occupation num¬ 
bers of single-particle states. Below we provide an exam¬ 
ple, in which the Lagrange multipliers associated with these 
conserved occupations enter as parameters in a generalization 
of the equilibrium transformation Eq. (3): more specifically, 
these multipliers play the role of effective inverse tempera¬ 
tures of the individual single-particle states. On the other 
hand, in non-integrable cases, the eigenstates of the Hamil¬ 
tonian are not single-particle states. Then one generically ex¬ 
pects the stationary state of the system to be in thermal equi¬ 
librium at a temperature T - IIP, which is determined by the 
initial conditions of the dynamics of the system. Accordingly, 
the enhanced symmetry that is present in the stationary state 
of the non-interacting integrable system breaks down and the 
corresponding Schwinger-Keldysh action is invariant under a 
single Tp, only for that specific value of p. This shows that the 
transformation Tp can be generalized in order to account for 
the appearance of a generalized Gibbs ensemble in the trivial 
case of a system that can be diagonalized in terms of single- 
particle states. However, the question whether the generalized 
Gibbs ensemble emerging in the stationary states of generic 
integrable systems is characterized by a symmetry involving 
the Lagrangian multipliers associated with the respective in¬ 
tegrals of motion as parameters, is beyond the scope of the 
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present work. 

As an example, let us consider bosons on a (7-dimensional 
lattice with nearest-neighbour hopping and on-site interaction 
(i.e., the Bose-Hubbard model [72]), with Hamiltonian 

^ ~ ^kin ^int» 

77kin — ^ ^ j 5 

(U'> (15) 

Hint = y X “l (“l “1 “ ’ 

where ai is the annihilation operator for bosons on the lat¬ 
tice site 1 , f is the hopping matrix element between site 1 
and its nearest-neighbours 1 ', while U determines the strength 
of on-site interactions. We hrst consider the case U - 0, 
which is trivially integrable: the kinetic energy contribution 
to the Hamiltonian is diagonal in momentum space and the 
corresponding single-particle eigenstates are the Bloch states. 
These are labelled by a quasi-momentum q, and in terms 
of creation and annihilation operators for particles in Bloch 
states, flq and Oq respectively, the kinetic energy can be writ¬ 
ten as 

^kin “ 'y ] ^q^q^q- (16) 

q 

Let us now consider a Schwinger-Keldysh functional integral 
description of the stationary state of the system. Then, the 
kinetic energy in Eq. (16) yields a contribution to the corre¬ 
sponding action which reads 

S «.kin = - (‘Aq,+iAq,+ " fq-^q-) , (17) 

'Jt q 

where i/(q_+ and tpq- are the helds on the forward and back¬ 
ward branches of the closed time path respectively, expressed 
in the basis of Bloch states. is invariant under the trans¬ 

formation of the helds 

Vq,^{uj) = (18) 

where, as in Eq. (3), we arrange the helds in a spinor 
'l'q,o-(w) = (iAq,o-(w), ^*q o-(-«)) . The cmcial point is that 
Pq can be chosen to depend on the quasi-momentum q, indi¬ 
cating that to each eigenstate of the system we can assign an 
individual “temperature” Tq - IjPq such that the correspond¬ 
ing mean occupation number Wq = {a\p.^ is determined by a 
Bose distribution with precisely this “temperature.” 

Let us now consider the opposite limit in which the hopping 
amplitude t vanishes while the interaction strength U is hnite. 
The interaction energy Hint in Eq- (15) is diagonal in the ba¬ 
sis of Wannier states localized at specihc lattice sites and the 
occupation numbers h\ - a^ai of these sites are conserved, 
rendering the system integrable. The generalized symmetry 
transformation appropriate for this case can be obtained from 
Eq. (18) by replacing the quasi-momentum q by the lattice site 
index 1 and by introducing a set of “local (inverse) tempera¬ 
tures” Ti (jS|) instead of Tq (Pq). 


In the generic case, when both the hopping t and the in¬ 
teraction U are non-zero, the system is not integrable. Then, 
neither the generalized transformation Eq. (18) nor its variant 
with local “temperatures” are symmetries of the correspond¬ 
ing Schwinger-Keldysh action, showing that this case eventu¬ 
ally admits only one single global temperature, which deter¬ 
mines the statistical weight of individual many-body eigen¬ 
states of the system. 


B. Dissipative contributions in equilibrium 

The functional integral with the action S in Eq. (8), as it 
stands, is not convergent but it can be made so by adding to 
S an inhnitesimally small imaginary (i.e., dissipative) contri¬ 
bution [55, 56]. Within a renormalization-group picture, this 
infinitesimal dissipation may be seen as the “initial value,” at 
a microscopic scale, of finite dissipative contributions, which 
are eventually obtained upon coarse graining the original ac¬ 
tion S and which result in, e.g., hnite lifetimes of excitations 
of the effective low-energy degrees of freedom. The precise 
form of the corresponding effective low-energy action and, 
in particular, of the dissipative contributions which appear 
therein, is strongly constrained by the requirement of invari¬ 
ance under Tp of the starting action at the microscopic scale; 
in fact, terms which violate this symmetry will not be gener¬ 
ated upon coarse-graining. In the discussion below we iden¬ 
tify those dissipative contributions to the Schwinger-Keldysh 
action which are invariant under 7]g. This allows us to an¬ 
ticipate the structure of any low-energy effective action pos¬ 
sessing this a symmetry. Note, however, thaf finite dissipative 
terms may appear even at the microscopic scale because of, 
e.g., the coupling of the system to an external bath. Below we 
consider two instances of this case; in Sec. VIB we show that 
Tp cannot be a symmetry of the action if the system is cou¬ 
pled to Markovian baths and driven — a situation described 
by a quantum master equation. Another specihc example, in 
which the equilibrium symmetry is realized, is the particle 
number non-conserving coupling of the Schwinger-Keldysh 
action Eq. (8) to an ohmic bath. This situation, which we dis¬ 
cuss in Sec. IV C, is of particular interest, because its classical 
limit renders what is known as the dynamical model A [63] 
with reversible mode couplings (termed model A* in Ref. 64); 
this correspondence allows us to establish a connection with 
the known equilibrium symmetry of the generating functional 
associated with this classical stochastic dynamics. 

Below we discuss dissipative terms of the action invariant 
under Tp, which involve hrst single particles (being quadratic 
in the helds of the Schwinger-Keldysh action) in Sec. IV B 1, 
and then their interactions in Sec. IV B 2. 


1. Single-particle sector 

Dissipative contributions to the single-particle sector of the 
Schwinger-Keldysh action which are invariant under Tp take 
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the form 

Sd^i I h(u, q) q)</)c{co, q) - (pg(aj, q)<^*(w, q) 

+2 coth(;Sw/2)0* (w, q)(l>q{aj, q)), (19) 

with an arbitrary real function h(a), q) which transforms under 
time reversal as h{ai, q) = /r(w, -q). When such dissipative 
terms are introduced in order to regularize the Schwinger- 
Keldysh functional integral, a typical choice for h(a),q) is 
h(a), q) - e [55, 56] with e —» 0. This ensures that the Green’s 
functions in the absence of interactions satisfy a fluctuation- 
dissipation relation (we postpone the detailed discussion of 
such relations to Sec. VIA). The FDR for non-interacting 
Green’s functions, together with the invariance of interactions 
under the transformation T^g shown in Sec. IV A 2, guaran¬ 
tees that the FDR is satisfied to all orders in perturbation the¬ 
ory [40]. 

While there are no restrictions on the form of the function 
h(ai, q), the hyperbolic cotangent coth(y6cu/2) appearing in the 
last term of 5;/ is uniquely fixed by the requirement of in¬ 
variance under Tp, as can be verified by following the line of 
argument presented in Appendix A. In particular, S d with a 
certain value of p in the argument of coth(y8w/2) is invariant 
under Tpi if and only if p' - p. This shows that, remarkably, 
the appearance of the thermodynamic equilibrium Bose distri¬ 
bution function niuS) - l/(e^" - 1) at a temperature T - Ijp 
in coth(/lw/2) = 2 n{ai) + 1, can be traced back to the fact that 
T^g is a symmetry of the action. Note that for simplicity we 
considered here only the case of vanishing chemical potential, 

- 0. For finite p, the frequency o) in the argument of the 
hyperbolic cotangent in Eq. (19) should be shifted according 
to tu —> w - /i, as we discussed in point 5 in Sec. III. 


2. Dissipative vertices 

The dissipative contributions discussed in the previous sec¬ 
tion are quadratic in the held operators and they naturally oc¬ 
cur, e.g., when the system is coupled to a thermal bath by 
means of an interaction which is linear in those fields. How¬ 
ever, this type of coupling necessarily breaks particle num¬ 
ber conservation. The number of particles is conserved if 
instead the system-bath interaction term commutes with the 
total number of particles of the system, N = n(x), where 

n(x) = i/t^(x)i/t(x) is the local density. In other words, to ensure 
particle number conservation, it is necessary that the coupling 
terms are at least quadratic in the system operators. Accord¬ 
ingly, dissipative vertices appear in the Schwinger-Keldysh 
action after integrating out the bath degrees of freedom. Then, 
the requirement of invariance of these terms under Tp allows 
us to infer a priori their possible structure. In particular, we 
And that a frequency-independent number-conserving quartic 
vertex (i.e., the dissipative counterpart to the two-body inter¬ 
action oc in the Hamiltonian (11)) is forbidden by the 
thermal symmetry. 

A generic quartic vertex, which conserves the number of 


particles and which is local in time, can be parameterized as 

Sd - —i I d(a»i — CU2 + W3 — 0 ) 4 ) 

*JCO] ,...,(04 

X [/i(wi,CU2,W3, 014 ) 1 /'* (wl)t/'+(W2)l/'+(W3)lA+(<^4) (20) 

+f 2 { 0 Ji,C 02 , 0»3, Ol4)l/'!(Wi)t/'-(W2)l/'* (W3 )iA-(<4^4) 

4-/3(011, W2, 0»3,014)1/'* (0li)t/'+(W2)l/'l(0l3)l/'_(0l4)], 

where /i_ 2,3 are real functions; in order to simplify the notation 
we do not indicate the (local) spatial dependence of the fields, 
which is understood together with the corresponding integra¬ 
tion in space. Conservation of particle number is ensured by 
the t/(l) invariance i/'± e“*i/'± on each contour separately, 

with generic phases a±, while the overall d-function on the 
frequencies guarantees locality in time. Restrictions on the 
functions /i 2,3 in the generic dissipative vertex in Eq. (20) 
follow from the requirements of causality [55], according to 
which Sd Co the action must vanish for 1 /'+ = ip-, and invari¬ 
ance of the dissipative vertex under the equilibrium transfor¬ 
mation. These conditions are studied in detail in Appendix B. 
In particular, we find that they cannot be satisfied if f\pp, are 
constant, i.e., do not depend on the frequencies. One particu¬ 
lar choice of these functions that is compatible with the con¬ 
straints is given by 

/l(wi, W2, W3, W4) = /2(wi,W2, W 3 ,W 4 ) 

= (wi - W2)coth(y3(wi - W 2 )/2), (21) 

/3(wi, W 2 , 0 ) 3 , 014 ) - -4 (oji - 0 J 2 ) (n(wi - 0 J 2 ) + 1), 

with the Bose distribution function n{a>). It is interest¬ 
ing to note that, in the basis of classical and quantum 
fields, this corresponds to a generalization of Eq. (19) with 
h(o>, q) = a>, in which the fields are replaced by the respec¬ 
tive densities defined as pc - ( 1 /'*/+ 4- /*/-) / V2 and pq - 
( 1 /'* 1 /'+ -(/'*(/'_)/ V 2 . Another notable property of this solu¬ 
tion is that for co\p —> 0 we have f\ 22 {o)i, 002 , 043 , 0 ) 4 ) 2T 
and f 3 (o>i,oJ 2 ,o> 3 ,o> 4 ) —» -4r, i.e., these limits of vanishing 
frequencies are finite. This implies that the form of Sd with 
/i, 2,3 given by Eq. (21) is to some extent universal: indeed, it 
should be expected to give the leading dissipative contribution 
to the Schwinger-Keldysh action of any number-conserving 
system in the low-frequency limit. At higher frequencies, 
other less universal solutions might also be important and one 
cannot make a general statement. 

C. Classical limit, detailed balance and microreversibility 

A transformation analogous to 7/ — which becomes a 
symmetry in equilibrium — was previously derived for the 
stochastic evolution of classical statistical systems in contact 
with an environment, within the response functional formal¬ 
ism [44-47, 73-76]. This formalism allows one to determine 
expectation values of relevant quantities as a functional in¬ 
tegral with a certain action known as response functional, 
which can also be derived from a suitable classical limit of 
the Schwinger-Keldysh action for quantum systems [55, 56]. 
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In these classical systems, the environment acts effectively as 
a source of stochastic noise over which the expectation values 
are taken. 

Here, we show that the classical limit of Tp [58] yields ex¬ 
actly the transformation which becomes a symmetry when the 
classical system is at equilibrium [48]. In order to consider 
this limit within the Schwinger-Keldysh formalism, it is con¬ 
venient to express the Schwinger-Keldysh action in Eq. (8) in 
terms of the classical and quantum fields (foe and (foq, respec¬ 
tively, dehned in Eq. (4), and to reinstate Planck’s constant 
according to [55, 56] 

S SIh, coth(y8m/2) — > coth(yS/i«/2), (foq ^ fi(foq. (22) 

Then, the action can be formally expanded in powers of ti in 
order to take the classical limit ^ > 0, and the classical part 

of the Schwinger-Keldysh action is given by the contribution 
which remains for H - 0. Note that the limit H ^ 0 consid¬ 
ered here is formally equivalent to approaching criticality in 
equilibrium at finite temperature T = >8 \ for which /3A —> 0, 
where A is the energy gap, which can be read off from the 
retarded Green’s function (see, e.g.. Ref. 77). This equiva¬ 
lence conforms with the expectation that quantum fluctuations 
generically play only a subdominant role in determining the 
critical behavior of quantum systems at finite temperature. In 
order to see the emergence of a stochastic dynamics driven by 
incoherent (thermal) noise from a quantum coherent dynam¬ 
ics, we supplement the Schwinger-Keldysh action in Eq. (8) 
(describing the latter) with dissipative terms arising from its 
coupling to a bath. Eor simplicity, we assume this bath to be 
characterized by an ohmic spectral density, while the system 
is assumed to have the Hamiltonian in Eq. (11). Deferring 
to Sec. VIB the discussion of the theoretical description of 
such a system-bath coupling, we anticipate here that the re¬ 
sulting contribution to the Schwinger-Keldysh action can be 
written as in Eq. (64), under the assumption that 7(w)v(w) is 
linear in the frequency, i.e., y(w)v(tu) = 2 koj and by choos¬ 
ing La-{co) ifocr{<JJ, q), with the thermal bath acting indepen¬ 
dently on each momentum mode [56]. Then, in the classical 
limit ^ > 0, we And 


S = 


I 
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V2 

(cTj. -t- iVl) idf H- 

2 m 


<1). -H i 2 KT^, 


■A f (tfof(foe<foq+C.C.'). (23) 

Jt,x ' 


This action has the form of the response functional of the equi¬ 
librium dynamical models considered in Ref. 63: it includes 
both a linear and a quadratic contribution in the quantum held 
(foq, but no higher-order terms. After having transformed the 
quadratic term into a linear one via the introduction of an aux¬ 
iliary held (which is eventually interpreted as a Gaussian ad¬ 
ditive noise), this quantum held can be integrated out and one 
is left with an effective constraint on the dynamics of the clas¬ 
sical held (foe, which takes the form of a Langevin equation; 
here: 


(i - k) d,(foe 




(foe + rj. 


(24) 


where rj = Tj(t,x) is a (complex) Gaussian stochastic noise 
with zero mean {r](t, x)) = 0 and correlations 

{T](t, x)Tj*(t',x')} = KT6(t - t')6^‘^\x - x'), (25) 

<;7(f,x)77(f',x')> = 0. (26) 

Equation (24) describes the dynamics of the non-conserved 
(complex scalar) held (foe without additional conserved densi¬ 
ties, which is known in the literature as model A [63]. How¬ 
ever, as can be seen from the complex prefactor i — k of the 
time derivative on the left-hand side of Eq. (24), the dynam¬ 
ics is not purely relaxational as in model A, but it has addi¬ 
tional coherent contributions, also known as reversible mode 
couplings [76]. The fact that the simultaneous appearance 
of dissipative and coherent dynamics can be described by a 
complex prefactor of the time derivative is specific to thermal 
equilibrium: in fact, dividing Eq. (24) by i - k, the reversible 
and iiTeversible parts of the resulting Langevin dynamics are 
not independent of each other and in fact their coupling con¬ 
stants share a common ratio [46, 47, 78]. Under more general 
non-equilibrium conditions, however, these reversible and ir¬ 
reversible generators of the dynamics have different micro¬ 
scopic origins and no common ratio generically exists. In the 
present equilibrium context, however, the action Eq. (23) cor¬ 
responds to model A* in the notion of Ref. 64, and the form 
of the classical transformation appropriate for this case which 
becomes a symmetry in equilibrium was given in Ref. 77. 
This transformation emerges as the classical limit of Tp dis¬ 
cussed in the previous sections [58]. In fact, foryS = —> 0 

and neglecting the contribution of the quantum fields in the 
transformation of the classical fields (i.e., at the leading order 
in Ti), Eq. (5) becomes 


7><l>c(f,x) = (Tx<A>e(-t,x), 

i i \ (27) 

Ti 3 %(t,x) = 0-;c(^<E>^(-f,x) -H —d,^e(-t,x)j, 

after a transformation back to the time and space domains. 
Upon identifying the classical held with the physical held 
and with the response held O, according to = ifl), 
Eq. (27) takes the form of the classical symmetry introduced 
in Ref. 48. Note, however, that the transformation (27) is not 
the only form in which the equilibrium symmetry in the clas¬ 
sical context can be expressed. In fact, the transformation of 
the response held O can also be expressed [46, 47] in terms of 
a functional derivative of the equilibrium distribution rather 
than of the time derivative of the classical held di<^e as in 
Eq. (27). The existence of these different but equivalent trans¬ 
formations might be related to the freedom in the definition 
of the response held, which is introduced in the theory as an 
auxiliary variable in order to enforce the dynamical constraint 
represented by the Langevin equation [46, 47, 76, 79] such 
as Eq. (24). This implies [79] that the related action acquires 
the so-called Slavnov-Taylor symmetry. As far as we know, 
the consequences of this symmetry have not been thoroughly 
investigated in the classical case and its role for quantum dy¬ 
namics surely represents an intriguing issue for future studies. 

We emphasize the fact that the derivation of the symmetry 
in the classical case involves explicitly the equilibrium prob- 
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ability density [46, 47]. Indeed, the response functional con¬ 
tains an additional contribution from the probability distribu¬ 
tion of the value of the fields at the initial time, after which 
the dynamics is considered. This term generically breaks 
the time-translational invariance of the theory [46, 47], un¬ 
less the initial probability distribution is the equilibrium one. 
Accordingly, when the classical equilibrium symmetry Tp is 
derived under the assumption of time-translational invariance, 
its expression involves also the equilibrium distribution. In 
the quantum case discussed in the previous sections, instead, 
time-translational symmetry was implicitly imposed by ex¬ 
tending the time integration in the action from -oo to -i-oo, 
which is equivalent to the explicit inclusion of the initial con¬ 
dition (in the form of an initial density matrix) and makes the 
analysis simpler, though with a less transparent interpretation 
from the physical standpoint. 

Although in classical systems this equilibrium symmetry 
takes (at least) two different but equivalent forms due to the 
arbitrariness in the definition of the response functional men¬ 
tioned above, it can always be traced back to the condition of 
detailed balance [46-48]. Within this context, detailed bal¬ 
ance is defined by the requirement that the probability of ob¬ 
serving a certain (stochastic) realization of the dynamics of the 
system equals the probability of observing the time-reversed 
realization, and therefore it encodes the notion of microre¬ 
versibility. This condition guarantees the existence and valid¬ 
ity of fluctuation-dissipation relations, which can be proved 
on the basis of this symmetry. In addition, detailed balance 
constrains the form that the response functional can take as 
well as the one of the equilibrium probability distribution for 
this stochastic process. 

The situation in the quantum case appears to be signifi¬ 
cantly less clear. In fact, a precise and shared notion of quan¬ 
tum detailed balance and quantum microreversibility is seem¬ 
ingly lacking. The first attempt to introduce a principle of 
quantum detailed balance dates back to Ref. 80, where it was 
derived from a condition of microreversibility in the context of 
Markovian quantum dynamics described by a Lindblad mas¬ 
ter equation. The mathematical properties of these conditions 
were subsequently studied in detail (see, e.g.. Refs. 81-85) 
and were shown to constrain the form of the Lindblad super¬ 
operator in order for it to admit a Gibbs-like stationary density 
matrix. However, even when this occurs, these operators are 
not able to reproduce the KMS condition and the fluctuation- 
dissipation relations because of the underlying Markovian ap¬ 
proximation, as we discuss in Sec. VIB. 

The notion of microreversibility in quantum systems ap¬ 
pears to have received even less attention, as well as its con¬ 
nection with some sort of reversibility expressed in terms of 
the probability of observing certain “trajectories” and their 
time-reversed ones. The definition proposed in Ref. 80 (also 
discussed in Ref. 86) appears to be a natural generalization 
of the notion in the classical case, as it relates the correla¬ 
tion of two operators evaluated at two different times with 
the correlation of the time-reversed ones. However, to our 
knowledge, the relationship between this condition and ther¬ 
modynamic equilibrium has never been fully elucidated. Al¬ 
though addressing these issues goes well beyond the scope of 


the present paper, they surely represent an interesting subject 
for future investigations. 


V. EQUIVALENCE BETWEEN THE SYMMETRY AND 
THE KMS CONDITION 

In this section we show that the invariance of the 
Schwinger-Keldysh action of a certain system under Tjg (as 
specified in Sec. IV) is equivalent to having multi-time corre¬ 
lation functions of the relevant fields which satisfy the KMS 
condition [42,43]. As the latter can be considered as the defin¬ 
ing property of thermodynamic equilibrium, this shows that 
the same applies to the invariance under the equilibrium sym¬ 
metry. 

The KMS condition involves both the Hamiltonian genera¬ 
tor of dynamics and the thermal nature of the density matrix 
which describes the stationary state of the system; heuris- 
tically this condition amounts to requiring that the many- 
body Hamiltonian which determines the (canonical) popula¬ 
tion of the various energy levels is the same as the one which 
rules the dynamics of the system. The equivalence proved 
here allows us to think of the problem from a different per¬ 
spective: taking the invariance under Tjj as the fundamen¬ 
tal property and observing that any time-independent Hamil¬ 
tonian respects it, we may require it to hold at any scale, 
beyond the microscopic one governed by reversible Hamil¬ 
tonian dynamics alone. In particular, upon coarse graining 
within a renormalization-group framework, only irreversible 
dissipative terms which comply with the symmetry (such as 
those discussed in Sec. IV B) can be generated in stationary 
state and the hierarchy of correlation functions respect ther¬ 
mal fluctuation-dissipation relations. The validity of KMS 
conditions (and therefore of the symmetry Tp) hinges on the 
whole system being prepared in a canonical density matrix 
p. Accordingly, if the system is described by a microcanon- 
ical ensemble, the KMS condition holds only in a subsystem 
of it, which is expected to be described by a canonical re¬ 
duced density matrix. Equivalently, this means that, in a mi- 
crocanonical ensemble, only suitable local observables satisfy 
this condition. In the case of quantum many-body systems 
evolving from a pure state, an additional restriction on the 
class of observables emerges due to the fact that, if thermal- 
ization occurs as conjectured by the eigenstate thermalization 
hypothesis (ETH) [8, 87, 88], the microcanonical ensemble is 
appropriate only if the observable involves the creation and 
annihilation of a small number of particles (low order cor¬ 
relation functions). This was shown to be also the case for 
EDRs [89, 90]: however, as pointed out above, the thermal 
symmetry implies the validity of EDRs involving an arbitrary 
number of particles, which leads to the conclusion that it does 
not apply to an isolated system thermalizing via the ETH. 
In other words, the thermal symmetry implies that the whole 
density matrix takes the form of a Gibbs ensemble, while in 
thermalization according to the ETH, only finite subsystems 
are thermalized by the coupling to the remainder of the sys¬ 
tem, which acts as a bath. Thus we see how Hamiltonian dy¬ 
namics favors thermal stationary states (with density matrix 
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p proportional to over arbitrary functionals p - piH). 
One explicit technical advantage of this perspective based on 
symmetry is that it allows us to utilize the toolbox of quantum 
field theory straightforwardly and to study the implications of 
Tp being a symmetry; this is exemplified here by considering 
the associated Ward-Takahashi identities and by showing the 
absence of this symmetry in dynamics described by Marko¬ 
vian quantum master equations in Sec. VIB. We also note that 
the presence of this symmetry provides a criterion for assess¬ 
ing the equilibrium nature of a certain dynamics by inspecting 
only the dynamic action functional, instead of the whole hi¬ 
erarchy of fluctuation-dissipation relations. In addition, this 
symmetry may be present in the actions of open systems with 
both reversible and dissipative terms. 

In the following, we consider a quantum system with uni¬ 
tary dynamics generated by the (time-independent) Hamil¬ 
tonian H, which is in thermal equilibrium at temperature 
T - and therefore has a density matrix p - j tre“^^. 
The KMS condition relies on the observation that for an op¬ 
erator in the Heisenberg representation A{t) - e'^'Ae one 
has 

A(t)p - pA(t — ip) (28) 

(for simplicity we do not include here a chemical potential, 
but at the end of the discussion we indicate how to account for 
it). This identity effectively corresponds, up to a translation of 
the time by an imaginary amount, to exchanging the order of 
the density matrix and of the operator A and therefore, when 
Eq. (28) is applied to a multi-time correlation function, it in¬ 
verts the time order of the involved times, which can be subse¬ 
quently restored by means of the quantum-mechanical time- 
reversal operation. Hence, the quantum-mechanical time re¬ 
versal naturally appears as an element of the equilibrium sym¬ 
metry Tp, while external fields have to be transformed accord¬ 
ingly, as indicated in Eq. (7). The application of time reversal 
yields a representation of the KMS condition which can be 
readily translated into the Schwinger-Keldysh formalism, as 
was noted in Refs. 38-40. In particular, it results in an infi¬ 
nite hierarchy of generalized multi-time quantum fluctuations- 
dissipation relations which include the usual EDR for two- 
time correlation and response functions of the bosonic fields 
as a special case (see Ref. 40 and Sec. VI A). One of the main 
points of this work is that these EDRs can also be regarded as 
the Ward-Takahashi identities associated with the invariance 
of the Schwinger-Keldysh action S under the discrete sym¬ 
metry [91] 7/? and that, conversely, the full hierarchy of EDRs 
implies the invariance of S under Tp. 

The argument outlined below, which shows the equivalence 
between the KMS condition and the thermal symmetry, in¬ 
volves several steps: as a preliminary we review in Secs. V A 
and V B how time-ordered and anti-time-ordered correlation 
functions can be expressed using the Schwinger-Keldysh tech¬ 
nique and we specify how these correlation functions trans¬ 
form under quantum mechanical time reversal. We apply 
these results to the KMS condition in Sec. V C: first we dis¬ 
cuss its generalization to multi-time correlation functions and 
then we translate such a generalization into the Schwinger- 
Keldysh formalism. This part proceeds mainly along the lines 


of Ref. 40, with some technical differences. Einally, we estab¬ 
lish the equivalence between the resulting hierarchy of EDRs 
and the thermal symmetry at the end of Sec. V C. 

A. Multi-time correlation functions in the Schwinger-Keldysh 
formalism 

a. Two-time correlation functions. Let us first consider 
a two-time correlation function 

{A{tA)B{tB)) = tr (A(tA)B(tB)p) (29) 

between two generic operators A and B (in the following we 
are particularly interested in considering the case in which A 
and B are the held operators i/r(x) or ij/^ (x) at positions x = x^ 
and X = Xb) evaluated at different times t^ and f^, respectively, 
in a quantum state described by the density matrix p. We as¬ 
sume that the dynamics of the system is unitary and generated 
by the Hamiltonian H. Then, the Heisenberg operator A at 
time f/i is related to the Schrodinger operator at a certain ini¬ 
tial time tj < tA via 

A(tA) = (30) 

with an analogous relation for B. 

The two-time correlation function can be represented 
within the Schwinger-Keldysh formalism as (see Appendix C) 

{A(tA)B(tB)} = {A.itA)B4tB)) 

= J lDmA_(tA)B4tB)e'^^'^K 

irrespective of the relative order of the times and tB- 
Here, the functional integral is taken over the fields 'T = 
, and the exponential weight with which a 
specific field configuration contributes to the integral is de¬ 
termined by the Schwinger-Keldysh action 5['i']. In the fol¬ 
lowing, by 0+/- we indicate that a certain operator O has 
been evaluated in terms of the fields i//± defined on the for¬ 
ward/backward branch of the temporal contour associated 
with the Schwinger-Keldysh formalism (see, e.g.. Refs. 55 
and 56). 

b. Multi-time correlation functions. We define multi¬ 
time correlation functions in terms of time-ordered and anti- 
time-ordered products of operators 

A(fA,l, ■ ■ .Ja.n) - a\{tA,l)ai{tA,2) ■ ■ ■ ClN(tA,N), 

B(tB,l, ■ Ib.m) - b M{tB,M)b M-\{tB,M-\) ' ' ' bi{tB,\), 

for f/ < tAA < ■ ■ ■ < tA,N < tf and f; < < ■ ■ ■ < 

tB,M < tf, where tf is an arbitrarily chosen largest time. Here, 
bosonic field operators. The specific sequence 
of time arguments in A and B (increasing and decreasing 
from left to right, respectively) leads to a time-ordering on 
the Schwinger-Keldysh contour: indeed, as we show in Ap¬ 
pendix C, the multi-time correlation function can be expressed 
as a Schwinger-Keldysh functional integral in the form 

(A(fA,l, . . . , tA,N)B(tB,\, ..., Ib.m)) 

- {B+(tB.l, ..., tB,M)A-{tA,\, . . .,tA,N))- (33) 
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c. Anti-time-ordered correlation functions. Not only 
time-ordered correlation functions such as Eq. (33) can be 
expressed in terms of functional integrals, but also correla¬ 
tion functions which are anti-time-ordered and which, e.g., 
are obtained by exchanging the positions of A(tA.\, ..., f^.w) 
and B(tB.i, ■.. ,tB,M) on the l.h.s. of Eq. (33). The construc¬ 
tion of the corresponding functional integral can be accom¬ 
plished with a few straightforward modifications to the pro¬ 
cedure summarized in Appendix C (and presented, e.g., in 
Refs. 55 and 56). In a stationary state one has [p,H] = 0 
and all the Heisenberg operators on the l.h.s. of Eq. (33) can 
be related to the Schrodinger operators at a later time tf. Then 
one finds 

{B(tB,l,. ■ ■ , tB,M)A{tA.\, ■ ■ ■ , fA.iv)) 

= (^+(fA.l> ■ • ■ j tA^N)B-{tBA, ■■■, tB,M))Sb^ (34) 

where the action S b describes the backward evolution and it is 
related to the action S which enters the forward evolution in 
Eq. (33) simply by a global change of sign St = S . 


B. Quantum-mechanical time reversal 

In this section we first recall some properties of the 
quantum-mechanical time reversal operation T [41] and then 
discuss its implementation within the Schwinger-Keldysh for¬ 
malism. T is an antiunitary operator, i.e., it is antilinear (such 
that — /!*T|i/') for A e C) and unitary (T^ = T“'). 

Scalar products transform under antiunitary transformations 
into their complex conjugates, i.e., (i/r|A|0) = {f/\A\f)*, where 
we denote by \f/) - Tli/r) and A - TAT* the state and the 
Schrodinger operator obtained from the state |i/f) and the op¬ 
erator A, respectively, after time reversal. Accordingly, ex¬ 
pressing the trace of an operator in a certain basis {\\p„)]n, one 
finds 

trA = = (trA)*. (35) 

n n 

In the last equality we used the fact that, due to the unitar- 
ity of T, also the time-reversed set {|i^„))„ forms a basis. Eor 
future convenience, we shall define the Heisenberg represen¬ 
tation of time-reversed operators such that it coincides with 
the Schrodinger one at time -tf, i.e., we set 

A(f,4) = (36) 

Note that this is distinct from the Heisenberg representation 
defined in Eq. (30), which coincides with the Schrodinger one 
only at time f,. In order to simplify the notation, we shall 
not distinguish these two different Heisenberg representations, 
assuming implicitly that the latter and the former are used, 
respectively, for operators and their time-reversed ones, such 
that A(f,) = A while A(-tf) - A. 

Let us now study the effect of time reversal on the generic 
multi-time correlation function in Eq. (33). Due to transla¬ 
tional invariance in time, the time arguments of the operators 


A and B can be shifted by f, - f/ without affecting the correla¬ 
tion function. Then, by using Eqs. (35) and (36), one has 

{A(tAf,. . . , tA.N)B{tB.\, ■ . . , tB.M)) 

- , -tA,N)B{-tB,l, ■ ■ ■ , -tB.M))}, (37) 

- {B\-tB,\, ■ ■ ■ ,-tB,M)A^ (-fA,l, • ■ ■ , -tA,N)}p, 

where the subscript in indicates that the expectation 

value is taken with respect to the time-reversed density op¬ 
erator p = TpT*, which is time-independent. The expectation 
value on the r.h.s. of Eq. (37) is anti-time ordered and there¬ 
fore it can be rewritten as a Schwinger-Keldysh functional in¬ 
tegral by using Eq. (34). The l.h.s., instead, is time-ordered 
and therefore it can be expressed as in Eq. (33), such that 
Eq. (37) becomes 

{B+{tBA, . . . , tB^M)A-{tA,l, ■ ■ ■ , tA,N)) 

- {A\(-tA,\,. ■ ■ , -tA,N)B*_{-tB,\, ■ ■ ■ , -tB,M))sb^ (3^) 

where the subscript bin St indicates that the sign of the action 
which describes the Hamiltonian evolution on the r.h.s. of this 
relation has been reversed, as explained below Eq. (34). The 
time-reversed action S differs from the action S associated 
with H which enters (implicitly, cf. Eq. (31)) Eq. (33) because 
in S the time evolution is generated by H, the initial state is 
the time-reversed density matrix p, and the integration over 
time extends from -tf to -f,. This latter difference becomes 
inconsequential as f,- ^ -oo and tf oo. 

Let us now consider the case in which A and B are products 
of the bosonic field operators i// and such that A± and B± 
involve the corresponding products of i//± and their complex 
conjugates. As there are no further restrictions on A and B, 
the l.h.s. of Eq. (38) can be genetically indicated as {(9['l']), 
where (9['l'] is the product of various fields on the Schwinger- 
Keldysh contour corresponding to B+{.. .)A_(...) which, ac¬ 
cording to the notation introduced in Sec. Ill, are collectively 
indicated by 'T = (i/r+, i/r*, (/r_, i/r* )^. With this shorthand no¬ 
tation, Eq. (38) can be cast in the form 

(Om) = <f?[T'l'])j^, (39) 

where the transformation 

T'T^(f,x) = 'Pl^(-f,x), (40) 

implements the quantum-mechanical time reversal within the 
Schwinger-Keldysh formalism. (With a slight abuse of no¬ 
tation, the same symbol T is used to indicate here both the 
quantum-mechanical time-reversal operator introduced above 
and the transformation of fields on the Schwinger-Keldysh 
contour in Eqs. (39) and (40).) In Eq. (40) we took into ac¬ 
count that the bosonic field operators in the Schrodinger pic¬ 
ture and in the real-space representation are time-reversal in¬ 
variant, i.e., = Ti/r(x)T* = i/^(x), which allows us to drop 

the tilde on the transformed field on the r.h.s. of Eq. (40). 
However, we note that in the last line of Eq. (37) the Her- 
mitean adjoint operators of those on the l.h.s. appear and this 
is the reason why both the r.h.s. of Eq. (38) and the trans¬ 
formation prescription Eq. (40) involve complex conjugation 
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of the fields. Note that the time-reversal transformation in 
Eq. (40) is actually linear, i.e., the complex conjugation af¬ 
fects only the fields variables and not possible complex pref¬ 
actors. This follows again from the last line of Eq. (37); the 
combination of two antilinear transfoimations (time reversal 
and Hermitean conjugation) results into a linear one. 

Let us mention that while Eq. (38) follows from the second 
line of Eq. (37), one could have equivalently taken its first line 
as the starting point for deriving a Schwinger-Keldysh time- 
reversal transformation. Then one would have been lead to a 
different implementation T' of time reversal: 

r'T<,(f,x) = 'T<,(-f,x), (41) 

with an additional overall complex conjugation of the corre¬ 
lation function. In some sense, the transformation T' is closer 
to the common way of representing the quantum-mechanical 
time reversal than T is, as it amounts to a mapping f i-> -f 
and i i-> -i [41]. Eor our purposes, however, T is of main in¬ 
terest, since it is part of the equilibrium transformation as we 
describe below. 

We finally note that an alternative implementation of the 
time-reversal transformation is based on the observation that 
the forward and backward branches of the closed-time path 
integral are actually equivalent if the dynamics of a system 
is time-reversal invariant (TRI) [55]. More specifically, this 
means that the Schwinger-Keldysh action is invariant — up to 
a global change of sign — upon exchanging the correspond¬ 
ing fields, i.e., This transformation 

is partly recovered in Eqs. (39) and (40): in fact, T in Eq. (40) 
involves an exchange of the contour indices, while the corre¬ 
sponding global change of sign in the action is indicated on 
the r.h.s. of Eq. (39) by the subscript Sb (cf. the definition of 
Sb below Eq. (34)). However, the time reversal transforma¬ 
tion T in Eq. (40) — derived from the quantum-mechanical 
time-reversal operation — additionally involves both complex 
conjugation and the time inversion f i-> -f. 

C. KMS condition and generalized fluctuation-dissipation 
relations 

The discussion of the previous sections about multi-time 
correlation functions and the time-reversal transformation 
provides the basis for the formulation of the KMS condition 
for multi-time coii'elation functions and of its representation 
in terms of Schwinger-Keldysh functional integrals. Eor the 
specific case of a four-time coiTelation function, the KMS 
condition is pictoiially illustrated in Eig. 1. As described in 
the caption, panel (a) summarizes the convention as far as 
the forward/backward contours are concerned. Panel (b), in¬ 
stead, refers to the KMS condition which, as anticipated af¬ 
ter Eq. (28), involves a contour exchange of the multi-time 
operators A and B. This exchange is the reason why the ar¬ 
rows in the second equality in Eig. 1 (b) are reversed, since 
both A and B turn out to be anti-time ordered when moved 
to the opposite contour. The appropriate time-ordering can 
be restored by means of the quantum-mechanical time rever¬ 
sal transformation introduced in Sec. V B, as indicated by the 


third equality in the figure. This is a crucial step, because 
only time-ordered correlation functions can be directly trans¬ 
lated into the functional integral by means of the usual Trot¬ 
ter decomposition, which makes the time-reversal transforma¬ 
tion indispensable in the construction. However, this does not 
mean that properties related to equilibrium conditions such as 
fluctuation-dissipation relations are fulfilled only if the Hamil¬ 
tonian is invariant under the time-reversal transformation. In¬ 
deed, it turns out that multi-time EDRs always involve both 
the Hamiltonian and its time-reversed counterpart [40], while 
as we show in Sec. VIA the EDR for two-time functions can 
be stated without reference to the time-reversed Hamiltonian, 
even if the Hamiltonian is not TRI. 

The KMS condition for a two-time function reads [92] 

{A{tA)B{tB)) = {BitB - i/3/2)A(tA + i/3/2)). (42) 

This relation can be proven by writing down explicitly the ex¬ 
pectation value on the l.h.s. with p - / tre~^^ and by in¬ 

serting the definition of the Heisenberg operators reported in 
Eq. (30). The generalization of this procedure to the case of 
multi-time correlation functions is straightforward and yields 

■ j tA,N)B(tB,l, ■ ■ ■ , Ib.m)} 

- {B(tB,i - 1/3/2 ,..., Ib.m - 1/3/2) 

X ^(fA,i + 1/3/2 ,..., tA.N + i/3/2)). (43) 

The real parts of the time variables on the r.h.s. of this equa¬ 
tion are such that the corresponding product of operators is 
anti-time-ordered (see Eig. 1). According to their definition in 
Eq. (32), A and B correspond to products of operators with, 
respectively, decreasing and increasing time arguments from 
right to left. Consequently, Eq. (43) can be expressed as a 
functional integral by using Eqs. (33) and (34) on the l.h.s. and 
r.h.s., respectively. The presence of an imaginary part in the 
time arguments of Eq. (43) does not constitute a problem; in 
fact, the functional integral along the vertical parts of the time 
path in Eig. 1 can be constructed by the same method as the 
horizontal parts, which is summarized in Sec. V A. Hence, we 
find 

{B+(tB.l, ■ ■ ■ , tB,M)A-(tA,l, ■ ■ ■ , tA,N)) 

— (^+(Ia,i + i/i/2 , ■.., Ia.n + i/i/2) 

X B-(tB.i - i/i/2,..., Ib.m - i/3/2))st. (44) 

As we did in Eq. (39) for the case of quantum-mechanical 
time-reversal, we may rewrite this equation in the form 

(Om) = {0['Kp'¥])s„ (45) 

where we define 

= T'_,(f - i(r/3/2). (46) 

This transformation 'Kp can be combined with the quantum 
mechanical time reversal T defined in Eq. (40) in order to ex¬ 
press the equilibrium transformation Tjj as 7]g = T o 'Kp. [93] 
By using Eq. (39) on the r.h.s. of Eq. (45), one concludes that 
the KMS condition implies 

{Om)^{0[Tp^])s, (47) 
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Figure 1 . (Color online) (a) Schematic representation of the Schwinger-Keldysh partition function [ 55 , 56 ]. The time evolution of the density 
matrix p(r) = can be represented by introducing two time lines to the left and right of p. These time lines correspond to the + and - 

parts of the Schwinger-Keldysh contour, respectively, (b) Schematic representation of the KMS condition for a four-time correlation function 
(«i(4,i)a2(4,2)62(4,2)61(4.1)) with 4j < 42 and 41 < 4,2, where <212 and 612 are bosonic field operators. As illustrated by the first equality 
(light blue box), this correlation function is properly time-ordered and therefore it can be directly represented within the Schwinger-Keldysh 
formalism with the operators ai,2 and bip evaluated along the - and + contours, respectively. The thermal density matrix p = can 

be first split into the products of x and then these two factors can be moved in opposite directions along the two time lines, with 

the effect of adding +ipi 2 and - 1 / 3/2 to the time arguments of ai,2 and bip, respectively. After these two factors have been moved to the end 
of the timelines, due to the cyclic property of the trace, they combine as represented by the second equality (orange box), where the time lines 
now take detours into the complex plane and the overall time order is effectively reversed as indicated by the arrows, which converge towards 
p instead of departing from it as in the case of sketch {a) or of the first equality of sketch ( 6 ). The original time ordering can be then restored 
by means of the time-reversal operation T, upon application of which operators are replaced by time-reversal transformed ones, p = TpT^ etc., 
and the signs of time variables are reversed. In addition, due to the anti-unitarity of T one has to take the Hermitian adjoint of the expression 
inside the trace. As a result, the order of operators is inverted and one obtains the third equality (green box) which is again properly time 
ordered. This construction can be generalized to arbitrary correlation functions, leading to Eq. ( 47 ). 


which indeed provides a generalized FDR for the correlation 
function ((9['f']) [40]. For various choices of the observable 
C)['F] we obtain the full hierarchy of multi-time FDRs, which 
contains as a special case the usual FDR for two-time func¬ 
tions (see Sec. VIA). Before demonstrating that this hierar¬ 
chy is actually equivalent to the invariance of the Schwinger- 
Keldysh action as expressed by Eq. (7), several remarks are in 
order; 

1. Although by means of the time-reversal transformation 
T we were able to restore the time ordering in Eq. (44), 
Eq. (47) still involves the time-reversed action S and 
not the original action S. However, in practice it will 
typically be clear how S can be obtained from S, e.g., 
by reversing the signs of external magnetic fields. In the 
absence of fields which break time-reversal invariance 
the tilde in Eq. (47) may be dropped, i.e., S - S. 

2. Equation (47) provides a generalized EDR expressed in 

terms of products (9['l'] of fields on the forward and 
backward branches, which we collected in the four- 
component vector T* = . A more fa¬ 


miliar formulation of EDRs is provided in terms of clas¬ 
sical and quantum fields (O^ and see Eq. (4)), which 
allow one to identify correlations functions (i.e., ex¬ 
pectation values involving only classical fields) and re¬ 
sponse functions or susceptibilities (expectation values 
involving both classical and quantum fields). EDRs pro¬ 
vide relations between correlation and response func¬ 
tions. In order to express the KMS condition in Eq. (47) 
in terms of the classical and quantum fields and 

(or, alternatively, of O = y>c, <Pq, <Pq) ), we note that 
they are linearly related to '!'+ and and therefore a 
generic product (9[<1>] of such fields can be expressed 
as a linear combination (with real coefficients) of prod¬ 
ucts (9, ['P]. According to Eq. (47), the expectation value 
of such a combination and therefore ((9[0]) can be ex¬ 
pressed as on its r.h.s. in terms of the same linear com¬ 
bination of since the transformation Tp is 

linear, it immediately follows that this linear combina¬ 
tion is nothing but {0{TpQ>])§, where the explicit form 
of the transformation of the components of the field <I) 
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under 7^ is reported in Eq. (5). The KMS condition 
then becomes 

= (48) 

In Sec. VIA, on the basis of Eq. (48), we derive the 
typical form of the EDR, which involves the correlation 
function of two classical fields and the susceptibility ex¬ 
pressed as a correlation between one quantum and one 
classical field. 

3. In the grand canonical ensemble with density matrix 
p = g-fi(fi-fiN) j ^ where V = 

is the particle number operator, the KMS condition 
Eq. (43) has to be generalized. In order to derive 
it, we split again the density matrix into a product 
^ jjjg caption of Eig. 1). 

Then, moving one of the two factors through each of 
the blocks of operators A and B as in the second (or¬ 
ange) box in panel (b) of Eig. 1 has not only the effect 
of adding +il3l2 and -ipjl to the time arguments of the 
field operators in A and B respectively, as in the case 
/i = 0; taking into account the canonical commutation 
relations, additional factors appears due to the fact that 

where cr - +1 and -1 for operators which are part of the 
time-ordered and anti-time-ordered blocks of operators 
A and B, respectively. The factors can be taken 

out of the expectation value in Eq. (43) (corresponding 
to the trace in Eig. 1 (b)) and do not affect the restoration 
of time order by means of the time-reversal transforma¬ 
tion, which is illustrated in the last (green) box in panel 
(b) of Eig. 1. Therefore, on the r.h.s. of Eq. (47) they 
would appear as prefactors, which are absorbed in the 
modified transformation given in Eq. (6). 


D. From the KMS condition to a symmetry of the 
Schwinger-Keldysh action 

In the previous section we showed that the KMS condition 
within the Schwinger-Keldysh functional integral formalism 
takes the form of Eq. (47) (with Tp given by either Eq. (3) or 
(6)). Here we argue further, that the latter relation is equiv¬ 
alent to requiring the invariance of the Schwinger-Keldysh 
action under the equilibrium symmetry Tp. To this end, we 
express the expectation values on the left and right hand sides 
of Eq. (47) as the functional integrals 

J 2)[T']<9['l']e'‘^''*'> (50) 

and 

(0[Tp'i'])g = J (51) 


respectively. Performing a change of integration variables 
T* —» 7p'¥ in the last expression, the argument of O simpli¬ 
fies because TpV —> = 'P, since Tp is involutive (see 

Sec. III). In addition, we show in Appendix D that the abso¬ 
lute value of the determinant of the Jacobian S' - 6{Tp'V)l6'V 
associated with Tp equals one, i.e., |Det,)7'| = 1, and there¬ 
fore the integration measure is not affected by the change of 
variable. Accordingly, one has 

{0[TpT])s = J (52) 

and by comparing this expression to Eq. (50) a sufficient con¬ 
dition for their equality is indeed Eq. (7), which expresses 
the invariance of the Schwinger-Keldysh action under the 
equilibrium transformation. Since the observable (9['T] in 
Eqs. (47), (50), and (52) is arbitrary, the condition is also 
necessary, which proves that Eqs. (7) and (47) (and, conse¬ 
quently, the KMS condition) are equivalent. 


VI. EXAMPLES 

In this section we discuss some concrete examples of how 
the invariance of a certain Schwinger-Keldysh action under 
the equilibrium transformation Tp can be used in practice. 
Eirst we show that Eq. (47) (or, equivalently, Eq. (48)) con¬ 
tains as a special case the quantum EDR which establishes 
a relationship between the two-time correlation function of 
the field (f>c (see Eq. (4)) and its response to an external per¬ 
turbation which couples linearly to it. This was also noted 
in Ref. 40; however, the conceptual advance done here con¬ 
sists in realizing that the EDR can be regarded as a Ward- 
Takahashi identity associated with the equilibrium symmetry. 
In Sec. VIB, instead, we elaborate on the non-equilibrium na¬ 
ture of Markovian dynamics described by a quantum master 
equation (of the Lindblad form), which is seen to violate ex¬ 
plicitly the equilibrium symmetry. 

The case of a system which is driven out of equilibrium 
by a coupling to different baths is considered in Sec. VIC. 
In particular, we discuss a single bosonic mode coupled to 
two baths at different temperatures and chemical potentials. 
Einally, in Sec. VID we briefly list a number of additional 
applications of the symmetry. Some of those have been put 
into practice in the context of classical dynamical systems, 
and could be generalized to the quantum case with the aid of 
the symmetry transformation discussed in the present paper. 


A. Fluctuation-dissipation relation for two-time functions 

Considering the invariance of the Schwinger-Keldysh ac¬ 
tion under Tp in Eq. (5) as the defining property of thermody¬ 
namic equilibrium, the generalized EDR in Eq. (48) (which, 
as discussed above, is nothing but the Ward-Takahashi iden¬ 
tity associated with the symmetry), emerges as a consequence 
of equilibrium conditions. Then, from the generalized EDR in 
Eq. (48), the EDR for two-time functions [55, 56] can indeed 
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be derived as a special case. The latter reads 

G^(w,q) = (G*(w,q) - G'^(w,q))cothOSw/2) 

= ! 2 Im G*(w, q) co\h(J3ul2), 

where the Keldysh, retarded, and advanced Green’s functions 
G^, G*, and G'^, respectively, are related to expectation values 
of classical and quantum fields via 

q) ( 27 r/^’ 6{cj - - q') = q)C(cu', q')), 

/G*(cu, q) (27r/^' 6{cj - " Q') = q)<^;(w', q')), 

iG'^iu, q) (27r/+‘ 6 {u) - - q') = 

(54) 

Here we are assuming translational invariance in both time 
and space, which is reflected in the appearance of frequency- 
and momentum-conserving d-functions in the previous ex¬ 
pressions. The FDR valid in classical systems [48, 49] can 
be recovered from the quantum FDR Eq. (53) by taking the 
classical limit as described in Sec. IV C. Contrary to what 
one might suspect at a first glance from the appearance of the 
time-reversed action S in the generalized FDR in Eq. (48), 
the derivation of the EDR for two-time functions we present 
below is valid irrespective of whether the action contains ex¬ 
ternal fields which break TRI or not. 

Let us consider the identity Eq. (48) for specific choices of 
the functional G[<h]. In particular, by taking O to be equal 
to (f>q(f>g, the expectation value on the l.h.s. of Eq. (48) has to 
vanish due to causality [55] and therefore 

0 = {<Pq(aj,q)4>;(co\q')) = -q)7>^;(a;',-q'))s ■ 

(55) 

Upon inserting the expression of the fields transformed ac¬ 
cording to Eq. (5), one readily finds the EDR 

Gf(w, q) = (Gf(w, q) - G|(w, q)) coth(/3w/2) (56) 

with the time-reversed action S. By setting, instead, O equal 
to the product of two classical fields, the l.h.s. of Eq. (48) 
renders the Keldysh Green’s function G^ in Eq. (54); by using 
the explicit form of Tp in Eq. (5) and the EDR derived above, 
the r.h.s. of that equation coincides with Gj and therefore one 
concludes that 

Gf(tu,q) = GV,-q), (57) 

which expresses the transformation behavior of the Keldysh 
Green’s function under time reversal of the Hamiltonian. 
Einally, by replacing G[<1>] in Eq. (48) with the product 
(^c(w, q)0J(w',q) of one classical and one quantum held, the 
l.h.s. renders by definition the retarded Green’s function G* 
while the r.h.s. can be worked out as explained above. Taking 
into account Eq. (56), one can eliminate the Keldysh Green’s 
function G^ appearing on the r.h.s. in favor of the retarded 
and advanced Green’s functions G* and G^, respectively, and 
eventually finds 

Gf(w,q) = GV,-q). (58) 


This relation, together with its complex conjugate (which 
relates the advanced Green’s functions calculated from the 
original and time-reversed Hamiltonians, respectively) and 
Eqs. (57) and (56), yields the EDR (53). 


B. Non-equilibrium nature of steady states of quantum master 
equations 

In classical statistical physics, the coupling of a system to 
a thermal bath and the resulting relaxation to thermodynamic 
equilibrium are commonly modelled in terms of Markovian 
stochastic processes, which can be described, e.g., by suit¬ 
able Langevin equations with Gaussian white noise [76]. The 
Markovian dynamics of a quantum system, on the other hand, 
is described by a quantum master equation in the Lindblad 
form [59, 60] (or by an equivalent Schwinger-Keldysh func¬ 
tional integral). Under specific conditions on the structure of 
the Lindblad operator [85, 94], the stationary state of this dy¬ 
namics is described by a thermal Gibbs distribution, such that 
all static properties (equal-time correlation functions) are in¬ 
distinguishable from those in thermodynamic equilibrium. In 
spite of this fact, however, dynamical signatures of thermody¬ 
namic equilibrium such as the KMS condition (see Sec. V C) 
and the EDR (see Sec. VIA) are violated [95, 96]. This vio¬ 
lation can be traced back to the fact that the Markovian and 
rotating wave approximations, which are done in deriving this 
quantum master equation cause an explicit breaking of the 
equilibrium symmetry as we show in this section — i.e., al¬ 
though the system is coupled to a bath in thermodynamic equi¬ 
librium, the system itself does not reach equilibrium. Physi¬ 
cally, this can be understood by noting that the microscopic 
dynamics underlying an approximate Markovian behavior in 
this case is indeed driven. A typical example in the context 
of quantum optics is an atom with two relevant energy lev¬ 
els separated by a level spacing ojq and subject to an external 
driving laser with frequency v detuned from resonance by an 
amount A = v - wq ^ wq. Transitions between the ground 
and the excited state are made possible only by the driving 
laser, and the energy scale which controls the validity of the 
Markov approximation in the dynamics of the two-level sys¬ 
tem is set by cjq. The excited state is assumed to be unstable 
and it can undergo spontaneous decay by emitting a photon to 
the radiation held, which acts as a reservoir. This illustrates 
the combined driven and dissipative nature of such quantum 
optical systems. 

In order to investigate in more detail the effect of the 
Markovian approximation typically done in the driven context 
on the validity of the equilibrium symmetry T^g, we consider 
a system with a certain action S, whose degrees of freedom 
are linearly coupled to those of a thermal bath. By integrating 
out the latter degrees of freedom, a dissipative contribution to 
the original action is generated. In order to simplify the dis¬ 
cussion, we assume that the bath consists of non-interacting 
harmonic oscillators labelled by an index p (with cr 

referring to the branch of the Schwinger-Keldysh contour), 
with proper frequency which are in thermodynamic equi¬ 
librium at a temperature T = !//?. The Schwinger-Keldysh 
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action of the bath is then given by 


to be linear in the bath variables and have a strength 


Sb 




(59) 


where the Green’s functions t') - {b/j^a-iOb*^ o-'(^^)) 

the oscillators of the bath are fixed by requiring it to be in 
equilibrium and therefore they read [55, 56] 


G/(t,t') ^ -i{n(aj,) + 

G;+(f, f) = 6(t - t')G/(t, f) + 6(t' - t)G;-(t, a 
G--(U t') = 0(t' - t)G-\t, t') + d(t - t')G;-(t, t'). 


Here n(oS) - Ij {e^‘^ - l) is the Bose distribution function and 
6 {t) denotes the Heaviside step function, which is defined as 


0(0 = 



t > 0, 
t < 0. 


(61) 


The coupling S sh between the system and the bath is assumed 


/i 

-r(ov(o--L-(o^;-(o). (62) 

Here, L±(t) are associated with the quantum jump or Lindblad 
operators, which we assume to be quasilocal polynomials of 
the system’s bosonic fields resulting from normally 

ordered operators in a second quantized description (e.g., the 
simplest choice would be L±(0 = iA±(6x)). In order to sim¬ 
plify the notation, we do not indicate here the spatial depen¬ 
dence of the fields (both of the system and of the bath), which 
is understood together with the corresponding integration over 
space. We assume the harmonic oscillators which constitute 
the bath to be spatially uncorrelated. The Schwinger-Keldysh 
functional integral with total action S +Sb+S sh involving both 
system and bath degrees of freedom is quadratic in the latter 
and, therefore, the bath can be integrated out. The resulting 
contribution is 


S' = 




Ah) y{co)v(ctj) 





X 


iGt^itj') Gt-(t,t')\( L4t')\ 

^G-(f,f') GZ-(t,t')j\-L^(nj’ 


(63) 


which eventually sums to 5. In deriving this action, we made 
the additional assumption that the bath modes form a 

dense continuum with a spectral density v(w) = 
centered around a frequency wq, with a bandwidth (see fur¬ 
ther below for its interpretation). Then, sums of the form 
Zi/j T/j ■ ■ ■ can be approximated as integrals over frequencies 

Aa> 'Y(to)v(cL>) ■ ■ ■, where yicj) describes the frequency 
distribution of the oscillator strengths. Inserting the explicit 
expressions (60) for the bath Green’s functions into Eq. (63), 
we obtain 


dw / r“ dto' 

S' - -i I —y(aj)v(<u) n(w)L* (w)L-((u) + {n((jj) -i- 1)L* (w)L+(w) - I —— {[0(w' - oj) {n{a>) + 1) 

Jwo-i? 2n \ J_oo 27t 


+0(-a)' + a>)n(a>)] L* (a;')T+(w') + \6{-a>' + oj) (n(cj) H- 1) -i- 6 (a)' - a>)n(a>)] L*_(a)')L^(a>')} , (64) 


where La-(aj) describe the quantum jump operators in the fre¬ 
quency space, and 6 (a)) - i'P^ + n6(co) (where V denotes 
the Cauchy principal value) is the Fourier transform of 6 (f) in 
Eq. (61). The terms involving the principal value contribute 
to the Hamiltonian part of the total Schwinger-Keldysh action 
of the system S + S'. Assuming that the jump operators are 
quasilocal polynomials of the bosonic field operators of the 


system, Lo-(aj) transform under Tj?/ as the field operators, i.e., 

Tf^’LAca) = e-‘^f^'^I^Ll(aj) and 7 >-L;(w) = e‘^f^'‘''I^L^(aj)- 

(65) 

inserting these expressions in Eq. (64) one finds that 
the contour-diagonal terms (i.e., those proportional to 
L'‘^(aS)La-(a)) with cr = +1, which include, in particular, the 
above-mentioned principal value terms) are invariant due to 
frequency conservation (cf. the discussion in Sec. IV A 2). On 
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the other hand, for the contribution which is off-diagonal in 
the Schwinger-Keldysh contour one finds 

^ ^r(o^)y(co) 

X [n(w)e^'"L+(w)L*(w) -i- (n(oj) + 1) e^^'"L^(w)L;(w)]. 

( 66 ) 

If the value of f}’ matches the inverse temperature fi - I IT of 
the bath modes, encoded in n(w), it is easy to see that these 
terms are invariant under 7^ because n(w)e^" = n{ai) H- 1. In 
summary, one concludes that 5' [7^'F] = 5'['I'] and being also 
the action S of the system in isolation invariant under 7^, the 
same holds for the total effective action S H- 5' of the system 
in contact with the thermal bath. 

In order to understand the effect of the Markovian approx¬ 
imation in the driven context on the invariance under 7^, let 
us now consider Eq. (64) after this approximation has been 
done. In particular, in order for these approximations to be 
valid, we assume that one can choose a “rotating” frame in 
which the evolution of the system is slow compared to the en¬ 
ergy scales 0)0 and § which characterize the bath. This is pos¬ 
sible if the system is driven by an external classical field such 
as a laser, so that the frequency of the drive bridges the gap 
between the natural time scales of the system and those of the 
bath. Then, all jump operators in Eq. (63) may be evaluated 
at the same time f, since the integral kernel in Eq. (63) (i.e., 
the product of bath Green’s functions , density of states 
v(w) and oscillator strength distribution yjw), integrated over 
the bath bandwidth) differs from zero only within a correla¬ 
tion time Tc ~ I/&, which is assumed to be much shorter than 
the timescale over which La-(t) evolves in the rotating frame. 
Additionally we assume that the spectral density v(w) of the 
states of the bath and the corresponding coupling strength 
^jy{o)) to the system do not vary appreciably within the rel¬ 
evant window 0)0 — "O' < 0 ) < 0)0 + "O', such that one can set 
y(w)v(w) ^ y{o)o)v(o)o). As before, the terms in Eq. (64) 
which are diagonal in the contour indices are invariant un¬ 
der Tj 3 (see Eq. (65)) also after the Markovian approximation; 
accordingly, we focus on the off-diagonal terms in Eq. (66), 
which become 


= -ij(o^o)y(o)o) 

X OO 1 

— [nLlio))L4o)) + (n + 1) L'Lio))L^(o))], (67) 

DC 27r 

where h = n(o)o) is the occupation number of the bath modes 
at frequency wq. This makes it clear that -& acts as a high- 
frequency cutoff, whose precise value, under Markovian con¬ 
ditions, does not affect the physics. Applying the transforma¬ 
tion 7/j to the fields one has 

= -iy(o)o)y(o)o) 

X oo 

^ [h/"L+(w)L*(w) + (h+ l)e-f‘‘^L4o))Ll(o))]. 

oo 

( 68 ) 


In order for to be invariant under Tjj, this expression 

should be equal to 5'g: in Eq. (67), which requires 

- h + I for all values of the frequencies o) within the 
relevant region o)o - § < o) < o)o + &■ Clearly, this is not 
possible and therefore the equilibrium symmetry is explicitly 
broken by the Markovian approximation in the driven context. 


C. System coupled to different baths 


A simple way to drive a system out of equilibrium is to 
bring it in contact with baths at different temperatures and 
chemical potentials. In this case, a net flux of energy and par¬ 
ticles is established across the system, preventing it from ther- 
malizing and, consequently, causing a violation of the symme¬ 
try Tj}. This scenario occurs, e.g., in the context of quantum 
electronics in a quantum dot connected by tunnel electrodes to 
two leads, between which a hnite voltage difference is main¬ 
tained. Here, for simplicity, we consider a minimal bosonic 
counterpart of this system, constituted by a single bosonic 
mode — described by the fields {i/^±, ) coupled to two baths 

of non-interacting harmonic oscillators kept at different tem¬ 
peratures and chemical potentials. We show explicitly that the 
non-equilibrium nature of this setup is accompanied by a vi¬ 
olation of the equilibrium symmetry. The generalization of 
this argument to a multi-mode system or to a larger number of 
baths is straightforward. 


We consider two baths of non-interacting harmonic oscil¬ 
lators bj^y, where the label i - 1,2 denotes the bath to which 
the operator belongs, while v denotes the corresponding mode 
with frequency w/y. Each bath is assumed to be in thermo¬ 
dynamic equilibrium with different inverse temperatures ySi, 
P 2 and different chemical potentials jUi, p 2 - The Schwinger- 
Keldysh action S hj of each bath takes the same form as in 
Eq. (59). However, in the present case, the distribution func¬ 
tions tiiio)) entering the bath Green’s functions depend on the 
chemical potentials jj., as n/(w) = ljjj 
previous section, the bath is assumed to be coupled linearly to 
the system variables L±(f) which are quasilocal polynomials 
of the bosonic helds The dynamics of the system 

and the baths is then controlled by a functional integral with 
the total Schwinger-Keldysh action S +Sb,i +Sh, 2 +Ssb, where 
S is to the action of the system, i.e., the single bosonic mode, 
and S sb is the system-bath coupling. As in the previous sec¬ 
tion, an effective dynamics for the system’s variables can be 
obtained by integrating out those of the bath. This yields an 
effective action 
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5' = 



A I f^oo j / , 

— 7(w)(m(w)L1(w)L-(w) + (n(w) + l)L*(w)L+(w) - I —— \[6(a)'- u){h{Lo) + 1) 

2n \ J_^ 2 ti \ 

+Q{-(J + w)n(w)] L\{io')L+(u)) + [6{-aj' + ai) (n(aj) + 1) + 0(a)' - w)n(w)] L* (w')i'-(w')| 


(69) 


The action S' is formally similar to the one in Eq. (64), 
with the spectral density ■y(a>)v(aj) replaced by the sum of 
the spectral densities of the baths J(a>) - Ji(a)) + J 2 (a>), 
where 7;(w) = y,(cu)v,(w) (with v,(w) and 7 /(w) defined as 
in Sec. VIB, see after Eq. (63)). Analogously, the distribution 
function n(a)) of the single bath we considered in Sec. VIB 
is replaced by the average of the distribution functions of the 
two baths i - 1,2 weighted by the relative spectral densities, 
i.e.. 


n(a>) - 


Ji(a>) 


Ji(a>) + J 2 (a>) 


ni{a>) 


J2(a>) 


J\(a)) + J 2 (a)) 


n2(u)). (70) 


Now, we consider how the effective action (69) transforms un¬ 
der the thermal symmetry Tp. Since here we are explicitly 
considering the presence of chemical potentials, we will use 
the generalization of the symmetry in Eq. (6). As dis¬ 
cussed in Sec. VIB, since are quasilocal polynomials 

of the bosonic fields of the system, they transform under 
as 




Accordingly, the products La-(to)L’^(oj) are invariant under 
the symmetry and therefore the contour-diagonal part of S', 
which contains such terms, is invariant. On the other hand, 
the part of S' which is off-diagonal in the Schwinger- 

Keldysh contour (i.e., the first two terms on the r.h.s. of Eq. 
(69)) is modified as: 


X +OO j 

CO 

X [n(w)e^(""'^)L+(w)r (w) 

-H (n(w) H- 1) e’^''^'^’L_(w)L;(w)]. (72) 

Comparing this expression with Eq. (69), one readily sees that 
the invariance of this term under requires h(a))e^^‘^^^^ - 
h(co) + 1, and therefore is not invariant under 

unless the two baths have the same temperature and chemi¬ 
cal potential, i.e., Pi = /32 - /S and jui = fi 2 - A*- In this 
case, one can easily verify from Eq. (70) that the average dis¬ 
tribution function h(a>) is just the Bose-Einstein distribution 
h(a)) = l/(e^t"-Ai) _ 1 ) as a consequence, n(cu)e^t"-A') - 
n(a)) + 1. 

In conclusion, when the system is driven out of equilibrium 
by a net flux of energy or particles, induced by a difference be¬ 
tween the temperatures or the chemical potentials of the baths, 
the total action of the system is no longer invariant under 7^^, 


D. Further applications 


Among the various possible applications of the symmetry 
Tja, we mention here: 

a. Symmetry-preserving approximations. Properties of 
interacting many-body systems can usually be obtained only 
by resorting to certain approximations. Then, while EDRs for 
correlation and response functions can be established exactly 
in the absence of interactions, in an approximate inclusion of 
the latter one has to make sure that the EDRs are not bro¬ 
ken. In other words, the approximation should conserve the 
thermal symmetry. This requirement for classical statistical 
systems has been implemented in the mode-coupling theory 
of the glass transition in Ref. 97. 

The field-theoretic formalism provides the natural frame¬ 
work for studying the behavior of systems at long wavelengths 
and low energies by employing renormalization-group meth¬ 
ods. In any of these methods, an effective description which 
is obtained by integrating out fast fluctuations must have the 
same symmetries as those present at microscopic scales. Eor 
example, in the case of the functional renormalization group 
(for reviews see Refs. 98-102), this is achieved by choosing 
an ansatz to approximate the scale-dependent effective action 
which incorporates these symmetries. In this context, the clas¬ 
sical limit of the thermal symmetry discussed here has been 
used in functional renormalization group studies of model 
A [103] and model C [104]. The quantum thermal symme¬ 
try, instead, is analogously preserved by the ansatz for the 
effective action chosen in Ref. [105] (in the form of EDT), 
where the scale-dependent crossover from quantum to clas¬ 
sical dynamics is studied. Alternatively, one can devise ap¬ 
proximation schemes which are compatible with the equiva¬ 
lent KMS conditions, as discussed in detail in Ref. 40. Note, 
however, that in this work the KMS condition in the form of 
Eq. (47) is imposed on the scale-dependent Green’s functions 
(supplemented by the corresponding condition on the vertex 
functions). On the other hand, the symmetry constraint can 
directly be applied to the effective action itself, which is the 
generating functional of vertex functions and contains infor¬ 
mation on all correlation functions. 

b. Fluctuation relations. Another concrete example of 
the usefulness of the thermal symmetry is provided by the 
derivation of transient fluctuation relations [106, 107] for 
time-dependent particle transport in Ref. [58]. There, the sym¬ 
metry is generalized in order to account for the presence of a 
time-dependent counting field which probes the current flow¬ 
ing through the system. This generalized symmetry yields a 
relation analogous to Eq. (47) (however, formulated in terms 
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of the generating functional for correlation and response func¬ 
tions), from which, e.g., a fluctuation relation for the probabil¬ 
ity distribution of work done on the system and the transmitted 
charge can be derived. 

VII. CONCLUSIONS 

We demonstrated here that the Schwinger-Keldysh action 
describing the dynamics of a generic quantum many-body 
system acquires a certain symmetry Tp if the evolution occurs 
in thermal equilibrium. To a certain extent, this symmetry 
was discussed in Ref. 58 in the specific context of fluctuation 
relations for particle transport. We traced the origin of this 
symmetry back to the Kubo-Martin-Schwinger (KMS) condi¬ 
tion which establishes a relationship between multi-time cor¬ 
relation functions in real and imaginary times of a system in 
canonical equilibrium at a certain temperature. Fluctuation- 
dissipation relations are then derived as the Ward-Takahashi 
identities associated with Tp. Remarkably, in the classical 
limit, this equilibrium symmetry reduces to the one known 
in classical stochastic systems, where it was derived from the 
assumption of detailed balance. By comparing with this clas¬ 
sical case, important questions on the nature of equilibrium in 
quantum systems arise. In particular, while microreversibility 
and detailed balance of the dynamics are deeply connected to 
the notion of equilibrium in classical stochastic systems, an 
analogous relationship for quantum systems does not clearly 
emerge and surely deserves further investigation. 

The equilibrium symmetry Tp is expected to play a cru¬ 
cial role in the study of thermalization in quantum systems, 
in particular when combined with a renormalization-group 
analysis. In fact, on the one hand, it provides a simple but 
powerful theoretical tool to assess whether a certain system 
is able to reproduce thermal equilibrium. This can indeed 
be accomplished by a direct inspection of the microscopic 
Schwinger-Keldysh action (or of the effective one generated 


after integrating out some degrees of freedom, e.g., along a 
renormalization-group flow) which describes the dynamics of 
the system, rather than checking, for instance, the validity of 
the fluctuation-dissipation relations among various correlation 
functions. On the other hand, the equilibrium symmetry might 
be useful also in order to investigate or characterize possible 
departures from equilibrium and, in this respect, it would be 
interesting to consider the case in which the system evolves 
in a generalized Gibbs ensemble [6, 8, 65-71]. Finally, while 
we focussed here on the case of bosons, the extension of our 
analysis to different statistics, for instance fermionic and spin 
systems, represents an interesting issue. 
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Appendix A: Invariance of quadratic dissipative contributions 

In order to show explicitly the invariance of 5 in Eq. (19) 
under the transformation Tp in Eq. (5), it is convenient to as¬ 
sume that h{uj, q) as a function of the frequency cu has definite 
parity and to consider separately the cases of odd and even 
functions, /!o(-w,q) = -/!o(a),q) and he(-LO,(\} - /!e(w, q), 
respectively. The generic case follows straightforwardly by 
linear combination. We then consider how the actions 



X (<l)c(w, q) -H coth(;8w/2)(I)^(w, q)) (A 1) 

— where O is introduced in Sec. IV A 1, see Eq. (12) — with 
a certain p transform under a transformation Tpi [see Eq. (5)] 
of the fields, with a generic parameter p'. One finds 


SATp.^] = ie J (sinh(^'w/2)<I)J(-ru, q) + cosh(/l'cu/2)Oj(-cu, q)) cr, 

X |^cosh(/l'w/2)<l)c(-cu, q) - sinh(/l'cu/2)<l)^(-w, q) coth(;6w/2) sinh(/3'w/2)Oc(-w, q) + cosh(/3'w/2)0^(-w, q))j . (A2) 

Note that the terms involving solely the classical held spinor cancel each other only if p' - p. Otherwise, terms oc Ojfhc 
remain, which lead to a violation of causality [55]. Eor P' - P instead, we obtain 

SATp<^] = -ie J (- ^mhipuHWciu, q) + cosh(^cu/2)(l>J(ru, q)) 

X (sir)h(Pto/2) - coth(;Scu/2) cosh(jSw/2)) <l>^(w, q). (A3) 


By means of the identity and after some straightforward algebraic manipulations one 

eventually finds that Sd[Tp<^] - 5rf[<l)], i.e., that Sd (with a 
certain P) is invariant under Tp. 


sinh X - coth x cosh x - - Ij sinh x. 


(A4) 
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Appendix B: Invariance of dissipative vertices 


Appendix C: Representation of correlation functions in the 
Schwinger-Keldysh formalism 


As pointed out in the main text, a first constraint that has to 
be imposed on the functions /i_ 2,3 appearing in the dissipative 
vertex in Eq. (20) follows from the requirement of causality of 
the Schwinger-Keldysh action. The latter must vanish when 
ip+ - \j/- [55], which implies the condition 

I -U)2 + U)i- W4)l/f* (wi)l/r+(W2)lA+(W3)lA+(t^4) 

X [/l(wi, W2, W3, W4) + f 2 {u)i,(x> 2 , W3, W4) 

+/ 3 (Wl> W3><4'4)] = 0. (Bl) 

Now let us consider Eq. (20) with the transformed fields 
Tp'i'cr, i-S-, SdiTpV]. Requiring it to be equal to 5^['T], we 
hnd that the following conditions should be fulfilled; 


f\{pj\, (jJ2, C1J3, C1J4) — /l(W2, CJi, OJ4, C1J3) — 0 , 

f2(a)i,a>2, W 3 , W 4 ) - wi, W 4 , CU 3 ) = 0, (B2) 

U> 2 , 0 J 4 ) ~ ^^^/3(W2, tUl , CU4, CU3) =0, 

where we used the conservation of frequencies implied by the 
d-function in Eq. (20) to simplify the exponent in the last 
line. Specifically, the necessary conditions are that the ex¬ 
pressions on the l.h.s. of these relations should vanish when 
integrated over frequencies after having been multiplied by 
the corresponding combinations of fields in Eq. (20) and by 
the (5-function on frequencies. The relations in Eq. (B2) are, 
however, sufficient conditions for the equality of S and 

Sdm- 

To begin with, we investigate the possible existence of a 
frequency-independent solution of Eqs. (Bl) and (B2) for 
/i, 2 , 3 ; these two equations then imply 

/i = -/2 = constant and - 0. (B3) 

However, this solution can be seen to lack physical relevance 
for the following reason: any physically sensible dissipative 
contribution to the Schwinger-Keldysh action compatible with 
the thermal symmetry can be considered as originating from 
integrating out a thermal bath which is appropriately coupled 
to the system. Anticipating the discussion of Sec. VIB, we 
note that such dissipative contributions always involve terms 
which are not diagonal in the contour indices (cf. Eq. (63)). S d 
with /i, 2,3 given by Eq. (B3), however, is not of this form. In 
fact, inserting Eq. (B3) in Eq. (20) yields a vertex that is equal 
to the two-body interaction in Eq. (11) apart from an overall 
factor of i, i.e., such a vertex would originate from an imagi¬ 
nary two-body coupling in a Hamiltonian. Clearly, this would 
violate hermiticity, rendering the Hamiltonian unphysical. 

While this demonstrates that — as anticipated in the main 
text — a frequency-independent number-conserving quartic 
vertex is not compatible with equilibrium conditions, solu¬ 
tions of Eqs. (Bl) and (B2) do exist with depending on fre¬ 
quency. One particular solution is given by Eq. (21) of the 
main text. 


c. Two-time correlation functions. In order to derive 
the representation of a two-time correlation function in the 
Schwinger-Keldysh formalism reported in Eq. (31), we in¬ 
sert the explicit expressions (30) for the Heisenberg operators 
Aitjf) and Bits) in the trace which defines the l.h.s. of Eq. (31) 
according to Eq. (29). Then, by introducing an additional and 
arbitrary time tj such that f,- < tA.B < tf, and by using the 
cyclic property of the trace one can write 

{A(tA)B(tB)) = tT(^e^‘^(‘f-‘‘‘)Be-‘^''‘^~'‘''pe'^''^-'‘^ . 

(Cl) 

The evolution of the density matrix is adjoint to the evolu¬ 
tion of Heisenberg operators, i.e., p(f) = e~‘^‘pe'^'. Thus, the 
operator acting from the left (right) on the 

density matrix p corresponds to the evolution in time from t' 
to t. In the correlation function (Cl), the time evolution from 
f; to tf on the left (right) of p is intercepted by the operator B at 
time tB (A at time tA). In order to convert the r.h.s. of Eq. (Cl) 
into a path integral, the standard procedure (see, e.g.. Refs. 55 
and 56) to be followed consists in writing the exponentials 
of the evolution operators as infinite products of infinitesimal 
and subsequent temporal evolutions (Trotter decomposition), 
in-between of which one can introduce completeness relations 
in terms of coherent states carrying the additional label “-I-” on 
the left of the density matrix, and a on its right. These co¬ 
herent states are eventually labeled by a temporal index on the 
forward (-I-) and backward (-) branches of the close-time path 
which characterizes the resulting action. Correspondingly, the 
operators on the left and on the right of the density matrix (B 
and A, respectively, in Eq. (Cl)) turn out to be evaluated on 
the fields (i.e., coherent states) which are dehned, respectively, 
on the forward and backward branches of the closed time path 
and this yields immediately the equality in Eq. (31), where 
the ordering of the matrix elements A_ and B+ on its r.h.s. is 
inconsequential. Eor the sake of completeness, we note that 
the expression as a Schwinger-Keldysh functional integral of 
a two-time function is not unique; in fact, it is straightforward 
to check that, by rearranging operators in Eq. (Cl), one can 
equivalently arrive at 


{A{tA)B(tB)) = 


|<A+(fA)B+(fB)) 

\{A^(tA)B^(tB)} 


for tA>tB, 
for tA<tB. 


(C2) 


However, as discussed below, the choice of Eq. (31) naturally 
lends itself to a generalization to multi-time correlation func¬ 
tions. 

d. Multi-time correlation functions. The functional inte¬ 
gral on the r.h.s. of Eq. (33) relation can be constructed from a 
straightforward generalization of Eq. (Cl): after a reshuffling 
of the operators such that A and B appear respectively on the 
left and right of the density matrix — as explained above — 
the temporal evolution can be artificially extended from f, to 
tf and it is intercepted on the l.h.s. of the density matrix by 
operators hi,..., at times ..., tB.M and on the r.h.s. by 
operators ai,.. ., 0 ^ at times i,..., tA.N- Again, the result¬ 
ing expression for the correlation function can be converted 
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directly into a path integral by a Trotter decomposition of the 
subsequent evolutions and by inserting completeness relations 
in terms of coherent states carrying the label “+” correspond¬ 
ing to the forward contour on the l.h.s. of the density matrix 
and the label for the backward contour on the r.h.s., which 
eventually leads to Eq. (33). 


The eigenvectors v, and eigenvalues T, of the frequency- 
dependent part, i.e., the solutions of the equation 

a>')vi(a)') = d,v,(w), (D3) 

2n 

are 


Appendix D: Jacobian of the equilibrium transformation 

In order to prove that |Det JTI = 1 it is convenient to calcu¬ 
late the Jacobian associated with Eq. (3) in frequency and 
momentum space, which reads 

= (27r)"'d("'>(q-Hq')/(w,w'), (Dl) 


vi(w) = (0,0,-^^"''^l)^, V 2 (co) = (-^^"''^l,0,0)^, 

V3(w) = l)^ , V4(w) = 1,0,0)^ , 

(D4) 

with Ai - A 2 - -1, and Aj - A 4 - 1, so that Det7 = 
A 1 A 2 A 2 A 4 = 1. As for the momentum-dependent part of the 
Jacobian matrix Eq. (Dl), we note that its eigenvectors can be 
constructed with any function /(q) by taking the even and odd 
combinations /(q) + /(-q); 


where 


J(to, co') - 2n6(a> - co') 


0 

0 

0 0 


0 


0 


0 

0 

0 

-poijl 


0 ^ 
0 

g^£j/2 

0 


. (D2) 



+ q') ifiq') ± f(-q')) = + (/(q) + /(-q)). 


(D5) 

Thus the eigenvalues of this part are +1, and hence the abso¬ 
lute value of the Jacobian matrix is |Det^| = 1. 
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